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ABSTRACT 

It is shown that a dynamo can operate in an Active Galactic Nuclei (AGN) accretion 
disk due to the Keplerian shear and due to the helical motions of expanding and twisting 
plumes of plasma heated by many star passages through the disk. Each plume rotates 
a fraction of the toroidal flux into poloidal flux, always in the same direction, through 
a finite angle, and proportional to its diameter. The predicted growth rate of poloidal 
magnetic flux, based upon two analytic approaches and numerical simulations, leads to 
a rapid exponentiation of a seed field, ~ 0.1 to ~ 0.01 per Keplerian period of the inner 
part of the disk. The initial value of the seed field may therefore be arbitrarily small 
yet reach, through dynamo gain, saturation very early in the disk history. Because of 
tidal disruption of stars close to the black hole, the maximum growth rate occurs at a 
radius of about 100 gravitational radii from the central object. The generated mean 
magnetic field, a quadrupole field, has predominantly even parity so that the radial 
component does not reverse sign across the midplane. The linear growth is predicted 
to be the same by each of the following three theoretical analyses: the flux conversion 
model, the mean field approach, and numerical modelin g. The common featur e is the 



conducting fluid flow, considered in companion Paper I (IPariev &: Colgate! 120061 ) where 
two coherent large scale flows occur naturally: the differential winding of Keplerian 
motion and differential rotation of expanding plumes. 
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Introduction 



The need for a magnetic dynamo to produce and amplify the immense magnetic fields observed 
external to galaxies and in clusters of galaxies has long been recognized. The theory of kinematic 
magnetic dynamos has had a long history and is a well developed subject by now. There are nu- 
merous mon ographs and r e yiew art i cles d e voted to the magnet i c dyn a mos in astrophysic s , som e 
of which are: IParkeJ (ll97flh:lMoffat,lj dl97fih:lstb3 (E97^:ICowUnel dl98lh : [Roberts fc Sowardl dl992h: 



Childress et al 



199m: Zeldovich. Ruzmaikin. fc Sokolofij <jl983h : IPriestj <|l982l ): lBussel (jl99lh 



Krause fc Radler 



(jl98Cll ): iBiskampi (|1993l ): iMestell (119991 ) . Hundreds of papers on magnetic dynamos are published 
each year. Three main astrophysical areas, in which dynamos are involved, are the generation 
of magnetic fields in the convective zones of planets and stars, in differentially rotating spiral 
galaxies, and in the accretion disks around compact objects. The possibility of production of mag- 
netic fields in the central parts of the bla ck hole accretion disks in AGN has been pointed out by 



(Gailitis et al 



2000 



Chakra barti. Rosner. fc Vainshteinl ()1994l ) and the need and possibility for a robust dynamo by 



Colgate fc Lil (1997 ) . Dyn amos have been also observed in the laboratory in the Riga experiment 



2001 



and in Karlsruhe experiment (jStieglitz fc Miillerl l200ll ) . although these 



flows only partially simulate astrophysical ones. The flow resulting in a dynamo is essentially three 
dimensional flow and often, especially under astrophysical circumstances, is a chaotic or turbulent 
flow. 

The shear in a rotating conducting fluid amplifies the magnetic field in the direction p erpendic- 
ular to the shear and facilitates the growth of the magnetic field. Originally, iParkerl (|1955l ) proposed 
to combine the effects of kinetic helicity of the small scale motions of the fluid with the differential 
rotation to generate large scale magnetic fields in the Sun. Here we consider just such a dynamo in 
its application to the differentially rotating flow in the accretion disk around Central Massive Black 
Holes (CMBH) in the centers of galaxies. The necessary and robust source of helicity is provided 
by the rising and expanding plumes of the gas heated by the star passages through the accretion 
disk. This property of t he rotation of expanding plumes in a rotating frame is discussed at length 
in the companion paper. iPariev fc Colgate! (|2006l ). "A Magnetic aw Dynamo in AGN Disks: I. The 
Hydrodynamics of Star-Disk Collisions and Keplerian Flow" , which is ref erred to as paper I. This 
natural and unique coherent flow is supported by experimental evidence (jBecklev et al.l 120031 ) and 
is fundamental to the origin of a robust dynamo in an AGN accretion disk. 

The magnetic dynamo in the disk is the essential part of the whole emerging picture of the 
formation and functioning of AGNs, closely related to the production of magnetic fields within 
galaxies, within clusters of galaxies, and the still greater energies and fluxes in the inter-galactic 



medium. Bla c k hole formation, Rossby w ave torquing o f the accretion disk ([Lovelace et al. 



Li et al 



2000 



2001bl : IColgate et al.ll2003l ). jet formation (|Li et al.ll2001al ) and magnetic field redis- 



1999 



tribution by reconnection and flux conver sion, and finally particle acceleration in the ra dio lobes 



and jets are the key parts of this scenario ([Colgate fc Li 



1999 



Colgate. Li fc Pariev!l200lh . Finally 



we note that if almost every galaxy contains a CMBH and that if a major fraction of the free en- 
ergy of its formation is converted into magnetic energy, then only a small fraction of this magnetic 
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energy, as seen in the giant radio lobes (jKronberg et al.ll200ll ). is sufficient to propose a possible 



feed back in structure formation and in galaxy formation. 

This work is arranged as follows: in section[2]we briefly overview the ingredients of the star-disk 
collisions dynamo with a brief review of the disk conditions and star disk collisions from Paper I. 
In section [3] we introduce the flux conversion dynamo analysis with a discussion of the necessary 
reconnection and turbulence driven resistivity. In section [J] the mean field theory is developed, in 
section [5] the dynamo equations and numerical method are developed, and in section [6] the results 
of numerical calculations are presented in support of all three approaches. Finally, we end with the 
conclusions in section [71 CGS units are used throughout the paper. 



2. The Ingredients of the Star-disk Collisions Dynamo 

A poloidal magnetic field can be of two types distinguished by the reflectional symmetry in 
the equatorial plane: quadrupole (or even) and dipole (or odd). Quadrupole field has the same 
sign of the radial component above and below the disk plane. The radial component of the dipole 
field changes sign under the reflection in the disk plane, it vanishes exactly at the disk plane. 
Rigourous definitions and properties of the odd and even fields are given in Appendix A. As is 
evident in Figure [JA, the quadrupole field has a large radial component, both within and external 
to the disk and furthermore maintains the same radial direction in both spaces. On the other 
hand the differential shear of a dipole field, symmetric about the midplane and therefore with zero 
radial component, results in no winding of the flux within the disk and therefore no toroidal gain. 
Various higher multipoles than the quadrupole have an alternating radial component as a function 
of radius and therefore a greater possibility of cancellation by reconnection. Differential winding 
of a symmetric poloidal field by the Keplerian flow results in a uniform toroidal field having the 
same direction over the disk height and within the disk, Figure [TB- An a deformation resulting in 
a large scale helicity, on the scale comparable to the radius of the disk, will transform toroidal field 
into poloidal field. This transformed field, or new poloidal flux must have the same polarity as the 
original poloidal flux. Then the closure of the dynamo cycle demands that this transformed flux 
be merged or reconnected with the original poloidal flux in order that it is augmented and hence 
produce gain. If this transformed flux alternates in direction (as would be the case for a dipole 
field across the thickness of a disk), then the merged flux will be averaged to near zero. Only in 
the case of the quadrupole field is there a possibility of a coherent addition to the original poloidal 
field when the a deformation, as produced by star-disk collisions, changes sign across the midplane 
and further rotates only tt/2 radians. One notes that star collisions in the opposite, axial, direction 
equally contribute to the quadrupole poloidal flux. The toroidal field produced by the shear of 
differential rotation from the quadrupole field, (Figure [T]B), has opposite directions far above the 
surface of the disk from that inside the disk. The opposite direction of the toroidal field above the 
disk is not shown in these drawings, because, in addition to the dynamo, it presumes the formation 
of a force-free oppositely directed helix in the conducting half space above and below the disk. We 
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have, however, predicted and calculated this force-free helix (iLi et al.ll2001al ). and furthermore, as 
mentioned above, we associate partial dissipation of its free energy with the visible structure of AGN 
jets. However, the magnitude of the quadrupole field in the region closer to the disk surface and to 
the midplane should be stronger (as computations actually prove). Therefore, the a deformation 
will primarily take the bottom portion of the quadrupole flux and convert it into radial flux above 
the disk plane directed in the same way as the upper portion of the quadrupole field, Figure [Tp. 
Therefore, in the accretion disk dynamo, plumes from star-disk collisions entrain and rotate toroidal 
flux by ~ 7r/2 radians, originating primarily from within the disk, Figure HP . Furthermore these 
plumes terminate close to or not far above the surface of the disk, and so produce negligible rotation 
of flux not so displaced from the disk. We then expect this rotated flux, before rotating a further 
7r radians and so before self cancellation, to reconnect as loops of poloidal flux, Figure Hp. These 
loops of flux now merge with the initial poloidal field, Figure [JA, thereby completing the cycle. To 
proceed with the dynamo problem we need to utilize the following results from Paper I: 

1. The distribution of stars in coordinate and velocity space in the central star cluster of an 
AGN. 

2. The velocity, and density of the plasma in the disk and in the corona of the disk. 

3. The hydrodynamics of the flow resulting from the passage of the star through the disk, the 
plumes. 

We review briefly the properties of the plumes produced by the star disk collisions as they 
relate to the dynamo. Then with these results we estimate the conductivity in order to develop a 
flux rotation theory of the dynamo. 



2.1. The Untwisting or Helicity Generation by the Plume 

Let us first introduce a term which is used frequently below. Because of high conductivity of 
the plasma considered in this paper, the magnetic field is close to a "frozen- in" state, when the 
magnetic field lines follow the motions of the plasma. Imagine now a closed contour attached to 
the particles of plasma with some magnetic flux passing through this contour. Let us also draw this 
contour such that it is close to being a plane contour. As the result of plasma motions, this contour 
can be rotated by some angle. If this rotation happens quickly enough, so that no substantial 
magnetic field crosses the contour due to diffusion, the magnetic flux passing through this closed 
contour remains almost unchanged. The component of the magnetic field normal to the plane of 
the contour and averaged over the surface of the closed contour should have rotated by the same 
angle as the contour. We will name this process "flux rotation". 

We describe the number density of central stellar cluster and the kinematics of the stellar orbits 
in Paper I. The most important result of this consideration is that the rate at which stars cross the 
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unit area of the disk surface peaks at a radial distance of about 100r 9 to 200r g from the CMBH, 
where r g = 2GM/c 2 = 3.0 • 10 13 M 8 cm = 9.5 • 10 _6 M 8 pc is the gravitational radius of the CMBH 
and M$ is the mass M of the CMBH expressed in units of 10 8 solar masses: M$ = M/10 8 M©. 
The rate of stars-disk collisions closer to the CMBH than 100r 5 is depleted because of the tidal 
destruction of stars in the gravitational field of CMBH and because of the grinding of the star 
orbits into the accretion disk plane. This grinding occurs because of the action of a drag, which 
every star experiences on its passage through the accretion disk. After many passages this drag 
causes an inclined Keplerian orbit of a star to become coplanar with the disk plane and this star 
becomes trapped inside the disk. 

The physics and dynamics of a star-disk collision is also considered in Paper I. Here we briefly 
summarize the results of Paper I for the convenience of the reader. A star collides with the disk 
at a typical velocity of 5 • 10 3 km/s to 10 4 km/s. The velocity of escape from the surface of a 
solar like star is 600km/s. This is one order of magnitude smaller than the velocity of the star 
moving through the gas in the accretion disk. Also, the sound speed in the accretion disk at a 
radial distance of ~ 200r 9 is ~ 50km/s (see Appendix A in Paper I for details and more accurate 
numbers). Therefore, the gravitational field of the star itself does not influence a highly supersonic 
flow of gas onto the star. In this regard, the physics of a star-disk collision is radically different from 
the physics of the classical accretion process on either the moving or the resting star. The classical 
theory of accre t ion o f interstellar gas with zero a ngula r momentum onto stars was d eveloped in 
Bondi & Hovk] dl944h : iBondi. Hovle Lvttletor] (|l947h : bondil (Il952h : iMcCreal (|l953h . Since the 
peculiar velocities of stars in the Galaxy are much less than 600km/s and the sound speed in the 
interstellar material is also much less than 600km/s, the gravitational potential of a star dominates 
the dynamics of the accretion flow in the near proximity of a star. The radius of the gravitational 
capture of the gas is much larger than the radius of the star. Captured gas falls almost radially down 
to the star surface. The presence of a small asymmetry or non-homogeneity of the surrounding gas 
causes nonzero angular momentum, which strongly influence the dynamics of the accretion flow 
below the gravitational capture radius. 

The term "collision" rather than "accretion" is much more appropriate for the description of 
the interaction of a passing star with the accretion disk. Because of the high velocity of the star, 
the cross section of the interaction of a star with the gas is equal to the geometric cross section of 
a star. The high ram pressure of the incoming stream with the density ~ 10~ 8 to ~ 1CP 10 gem -3 

10 6 K and 



strips away the outer layer of a star. The underlying layers with the temperature 
density ~ 10~ 5 gcm~ 3 are exposed. This picture is completely different from the physics of the 
mixing of the radial acc retion s tream with the stellar (solar) atmosphere in the classical accretion 
theory as described by iHovld (|1949l ). A radiation shock is formed in front of the star and the 
channel of the hot gas is left behind the star. This channel expands sideways inside the accretion 
disk and heats the surrounding gas. The hot gas is subject to buoyancy force acting away from the 
equatorial plane of the disk. As a result of this force, two plumes rising from the two sides of the 
accretion disk are formed at the location of the star-disk crossing. Note, that the amount of gas in 
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the rising plumes and the size of the plumes are much larger than the initial mass and size of the 
hot channel made by the star. 

As explained in Paper I, the plume should expand to several times its original radius by the 
time it reaches the height of the order 2H above the disk surface, where H is the semi-thickness of 
the disk. The corresponding increase in the moment of inertia of the plume and the conservation 
of the angular momentum of the plume causes the plume to rotate slower relative to the inertial 
frame. From the viewpoint of the observer in the frame corotating with the Keplerian flow at the 
radius of the disk at the location of the plume, this means that the plume rotates in the direction 
opposite to the Keplerian rotation with an angular velocity equal to some fraction of the local 
Keplerian angular velocity depending upon the radial expansion ratio. Since the expansion of the 
plume will not be infinite in the rise and fall time of ir radians of Keplerian rotation of the disk, 
we expect that the average of the plume rotation will be correspondingly less, or A<p < tt or ~ ir/2 
radians. Any force or frictional drag that resists this rotation will be countered by the Coriolis 
force. Finally we note that kinetic helicity is proportional to 



ft = y(Vxv). (1) 

For the dynamo one requires one additional dynamic property of the plumes. This is, that 
the total rotation angle must be finite and preferably ~ n/2 radians, otherwise a larger angle or 
after many turns the vector of the entrained magnetic field would average to a small value and 
consequently the dynamo growth rate would be correspondingly small. This property of finite 
rotation, A<f> ~ tt/2 radians, is a natural property of plumes produced above a Keplerian disk. 

Thus we have derived the approximate properties of an accretion disk around a massive black 
hole: the high probability of star-disk collisions, the three necessary properties of the resulting 
plumes all necessary for a robust dynamo. What is missing from this description is the necessary 
electrical properties of the medium. 



3. The Flux Rotation Dynamo 
3.1. The Conductivity of the Disk and the Corona 



The dynamos producing large scale magnetic fields require a compromise between high and low 
conductivity. A poor conductor or an insulator will not allow the field to be dragged with the motion 
of the medium. Ohmic dissipation will cause the magnetic field to decay, and if sufficiently rapid, 
no dynamo will be possible. In the limit of very high conductivity, kinematic exponential growth of 
magnetic field has been predict ed to occur in the presence of random o r chaotic three dimensiona l 



motions of the medium (e.g., IZeldovich. Ruzmaikin. Sokolofil 11983 : iRoberts &; Sowardl Il992l ) 



The problem with analytic three dimensional motions is that being described as kinematic, they 
are reversible in the sense that little or no entropy is generated by the motions themselves. The field 
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can be unwrapped by a "non-Maxwell" demon following the line of force. Since the "demon" does 
not have to "throw away" any information in following the reverse path, no entropy is generated. 
The chaotic behavior in time of kinematic mathematically reversible motions does not create entropy 
because of time reversal invariance of the equations. The plumes from star disk collisions indeed 
occur randomly in time, but since the initial state is as random as the final state, no entropy is 
generated just due to the randomness in time of the plumes themselves. By comparison, if the 
initial state were a large amplitude coherent wave, then phase scrambling would indeed alter the 
entropy, but by a relatively small amount compared to a scattering process that leads to a Maxwell 
distribution. This lack of a change in entropy is then equivalent to laminar flow (without molecular 
diffusion) where mixing is reversible. By contrast, the randomness created by fluid turbulence is 
irreversible, satisfying a principle of maximizing the dissipation of the free energy of shear flow in a 
fluid. The plumes, although random in time, result in a coherent addition of poloidal flux, because 
every plume translates axially, expands radially, and rotates through nearly the same angle, ~ 7r/2, 
for every plume. 

The negative effect of turbulence on dynamo g ain has been docume nted in thre e major liq 



uid s odium dynamo experiments: Lyon, Cadarache ( Bourgo in et al.ll2002i ). Maryland (jSisan et al 



20041 ) . and Madison (j Nornberg et al.ll2006l ; ISpence et al.ll2006l ). all using the similar flow configura- 



tions . Numerous theoretical simulations of these flows, the Dudley- James flow in a sphere (jDudlev &; James 
1989) or similar von Karman flow in a cylinder (i.e., two counter-rotating radially converging and 

axially diverging flow s ) in the kinematic or l aminar limit have been perform ed: (Lyon, Cadarache) 

Bourgoin etaD (12002 'iJPetrelis et all (boO ^Mane et al.N200.ll) . (Maryland nPefflev. Cawthorne fc Lathrop 



(20o"ofh"lSw(rl (-1 a! 



(2001 



), and fMadison) iBavliss et al.1 (|200fih : KTConnell et all iboOfJ ). Thev all 
predict exponential dynamo gain at a critical magnetic Reynolds numbers, Rm cr ,t ~ 50. Yet the ex- 
periments give a null result, i.e., no exponential dynamo gain for experimental flows where Rm, 



exp 



> 



130 ~ 2.5Rm cr jj are achieved in the ex periments. These null results are interpre t ed as due to the 



negative effects of turb ulent diffusion (jBourgoin et al.l 12004 ; ISpence et al.l 120061 : I Nornberg et al 



2006: lLaval et al.ll2006T ). Our generalized interpretation of these result s is that turbulence behaves 
as an enhanced d iffusion of magnetic flux or an enhanced resistivity (jBoldvrev Cattaneoll2004l ; 
Pontv et al.ll2005l ). In these experiments the turbu lent velocity v tur b ~ 0.4 < v > where < v > 



is the average shear velocity (j Nornberg et al.l 120061 ) . Then the turbulence le ads to a decreased 



(2) 



conductivity or an enhanced resistivity as described by lKrause h Radlerl (jl980l ) as: 

- a ° 

aturb ~ l + 47T/3<7 /c 2 ' 

where o"o = c 2 /(4irrjo) is the conductivity of the fluid, rjo is the magnetic diffusivity of the same 
fluid, and aturb is the effective conductivity in the presence of turbulence. The constant (5 is derived 
from mean-field electrodynamics assuming isotropic turbulence: 



P ^ (w/3)fLi, 



(3) 



where r rorr is the mean correlation time of a turbulent fluctuation. Since the correlation time is 



an eddy turn over time, then r c , 



Lcorr/vturb-, where L corr is an eddy size. We then identify 
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LcorrVturb/3 = /? as a turbulent diffusion coefficient and the turbulent conductivity becomes the 
original conductivity decreased by the factor 1 + 47t/3cjo/c 2 . In the limit of a large /3 » rjo, the 
effective magnetic diffusivity then becomes just the turbulent diffusivity, r) e ff = (3 ~ L corr vt ur b/3- 
It is the combination of turbulent diffusivity with fluid resistivity or equivalently, effective restivity 
that determines dynamo gain. In unconstrained flows L corr becomes the dimension of the largest 
eddy that can "fit" in the flow or L corr = (d(ln < v >)/(ix) _1 ~ L/2 where L is the dimension of 
the shear flow. This larger effective resistivity then results in a smaller effective magnetic Reynolds 
number that determines dynamo exponential gain, Rm e jj = L < v > /r/ e // for a given < v > and 
L. Since v tur t, ~< v > /2, then rj e ff ~ (l/3)(l/4) < v > L and Rm e fj ~ 12. This value of Rm e jj is 
significantly smaller than the predicted value of Rm cr j£ ~ 50 for the Dudley-James or von Karman 
flows used in the current major experi ments. It is even smaller than Rm rr ^ ~ 17 for Ponomarenko 



flow used in Riga dynamo experiment (jPonomarenko 



1973 



Gailitis fc Freiberdll976l ). 



The critical threshold for dynamo gain, Rm cr ^ is determined from kinematic dynamo calcu- 
lations without enhanced turbulent resistivity. For bounded sheared flows, namely except for all 
but these special flows mentioned above, Rm cr j£ ~ 100. The question is whether there can be any 
exponential dynamo gain in any unconstrained shear flows. 

By way of confirmation, in the two dyna mo experiments tha t have demonstrated positive 



exponenti al dynamo gain, the Rig a experiment (j Gailitis et al 



periment (jStieglitz &: Miiller 



2001 



2000 



2001 



and the Karlsruhe ex- 



U) , turbulence was greatly constrained by the presence of a ridged 
wall(s) separating the counter flowing shear flows. It was therefore well recognized that these ex- 
periments did not represent astrophysical dynamos, but, on the other hand, strongly confirmed 
dynamo theory. We were therefore convinced that a natural constraint of turbulence must exist for 
the dynamos of astrophysics. 

There are at least five constraints of turbulence in shear flow that may alter the magnitude of 
turbulence as well as its isotropy. We list five of these constraints, expecting others to be identified: 

1. viscosity. 

2. a ridged wall. 

3. a positive outward gradient of angular momentum. 

4. a gradient of entropy in a gravitational field (e.g., the base of the convection zone of stars). 

5. delay in the onset of fully developed turbulence in unconstrained shear flow. 

Viscosity may inhibit all turbulence so that in this limit the properties of turbulence are not 
relevant. 



A ridged wall affects the magnitude of turbulence and its isotropy (the law of the walls) . 
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The gradients of angular momentum and entropy apply to astrophysical circumstances, where 
depending upon the presence of other instabilities, e.g. the magnetorotational instability or Rossby 
vortex instability, turbulence may be a small fraction of the average shear flow. 

Time dependence is similar to viscosity, in the case where the initiation time of the shear flow 
may be very short compared to the development time of the turbulence as in the case of the plumes 
driven by star-disk collisions. This limit leads to negligible levels of turbulence compared to the 
shear flow. 

We chose a gradient of angular momentum and time dependence of the shear flow, e.g., plume 
flow, as the probable mechanisms of constraint of turbulence for the most likely circumstances for 
producing an astrophysical dynamo; angular momentum as the circumstance for the accretion in 
massive black hole formation, and time dependence for the constraint of the transient period of the 
rise and fall of plumes in astrophysical circumstances. 

It is with this uncertainty of the role of turbulence in dissipating the magnetic flux as opposed 
to amplifying it that the current work was undertaken. Hence, when we found the possibility of 
a combination of (a), a near laminar shear flow, Keplerian flow, and (b), a repeatable, transient, 
non-turbulent source of helicity, could the possibility of a robust astrophysical dynamo become 
evident. Although the star-disk collisions are random in time, the flow, to first order is repeatable 
and therefore not turbulent. The flow resulting from a superposition of many plumes may be 
chaotic in time, but the superposition of many plumes, all with the same rotation, leads to a net 
rotated flux in the same direction. On the other hand, the vortices in anisotropic turbulence make 
an arbitrary number of turns and so the instantaneous mean value of rotated flux is proportional 
to the square root of the number of vortices. Thus we characterize the plumes as semi-coherent 
rather than a truly chaotic phenomena in which the entropy would be increased. In addition we 
discuss next the analysis in which turbulence may augment or possibly limit the "fast dynamo". 



The dyna mos with non- vanishing 



growt h rate in the limit of very high conductivity are called 



fast dynamos (jVainshtein Z eldovich 



is stretch-twist-fold process (jSakharov 



1972). A classical picture of the fa st dynamo mechanism 



19821 ; IVainshtein &: Zeldovich! 119720. There are strong in 



dications that t he fast dynamo action is typical for chaotic flows (|Lau &: Finnlll993l ; 



Finn et al. 1991 



Finn 



1992; 



However, in the kinematic stage of the dynamo, a sharp exponential decrease in 
some spatial scales of the magnetic field occurs. The magnetic field becomes concentrated in the 
narrow sheets or narrow filaments until the frozen-in picture becomes invalid for any conductivity. 
In the kinematic limit the thickness of these structures of strong magnetic field is estimated as 
51 ~ LRm -1 / 2 , where the magnetic Reynolds number, Rm = vL/r), v is the velocity of the con- 
ducting fluid, L is the characteristic dimension of the fluid, and rj is magnetic diffusivity due to 
finite resistivity. The growth of small scale fields invalidates the kinematic approximation begin- 
ning with the resistive scale and up to the larger scales. The structure and the spectrum of this, 
hydromagnetic, turbulent dynamo is expected to be different from the hydrodynamic turbulence 
because of the action of the magnetic forces. The size 51 of the smallest magnetic structures depends 
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on the details and properties of the hydr qmagnetic regime of the turbulence, which are still the 



subject of active deb ate in the literature (|lroshnikovl Il963l ; iKraichnanl Il965l ; iGoldreich Sridhar 



19951 ; iBoldvrevi 12006) . but is always much smaller than the large scale L by some positive power 
of Rm. In the limit of infinite conductivity, no flux can merge in an infinite time and hence, there 
can be no multiplication of flux at a scale of the system (large scale). The motions with nonzero 
helicity h at a large scale and the ability to reconnect is required to obtain the growth of the large 
scale fields and magnetic flux comparable to the growth rate of the small scale field. When the 
large scale field growth is at a rate comparable to the growth rate of the small scale field, the 
characteristic growth time is of order of the diffusion time = L 2 /r]. 

Yet our ionized disk of thickness H = 2.6 x 10 cm, velocity v K « 10 9 cms and resistivity 
7] ~ 10 7 cm 2 s _1 at 1 eV ~ 10 4 K temperature, results in Rm ~ 10 15 . This is a number so large as 
to preclude useful growth of the magnetic flux and the large scale fields in a Hubble time, which is 
much shorter than the diffusion time of the magnetic field i<K// = H 2 /rj = 10 20 s. 

Only by invoking the phenomena of turbulent resistivity, (above) can the existence of an 
accretion disk dynamo producing large scale magnetic fields be made convincing. Turbulent resis- 
tivity within the disk is likely to be due to the same turbulence that creates the a-viscosity of the 
Shakura-Sunyaev disk or the Rossby vortices of the RVI disk. Within the disk we expect turbulent 
diffusion of the magnetic flux to be the same as that of angular momentum, and thus proportional 
to the Shakura-Sunyaev parameter a ss (see paper I). 

Reconnection may be occurring within the turbulence leading to more rapid dissipation of the 
magnetic flux than the turbulent cascade alone. However, the force-free fields above (and below) 
the disk that are produced by the winding of the dynamo-produced large scale fields need not be 
dissipated until they are projected large distances away from the disk. 

Recognizing this lack of fundamental understanding, but that laboratory and astrophysical 
observations lead to the same order of magnitude for reconnection, we proceed with the assumption 
that a value of Rm ~ 200 approximates the magnetic diffusion within the disk. In what follows 
this parameter could be several orders of magnitude larger, but not much smaller and still result 
in an effective accretion disk dynamo. 



3.2. Estimates of the Dynamo Growth Rate: The Flux Rotation Dynamo 

With these values of plume size, frequency, and magnetic diffusivity as well as Keplerian flow, 
let us make some estimates of the threshold parameters and the growth rate of the auo dynamo, 
which has been outlined in a previous section. The approach developed below, in the rest of 
section El we call the flux rotation dynamo as opposed to the mean field dynamo. Later we will 
compare these approaches with emphasis upon the difference between coherent motions versus 
random averaged variables. We consider for now the linear growth, i.e. when the magnetic field is 
not strong enough, such that one can neglect the back reaction of the generated magnetic field on 
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both the Keplerian flow and on the plume flow fields. 

Suppose that at some moment of time we have an even symmetry poloidal magnetic field, Bp 
(see Appendix [X] for definitions of even and odd symmetries). The radial component of this field 
within the disk defines a poloidal flux, Fp, such that at a given radius r, the poloidal flux through 
one half of the disk, either side of the mid-plane, is Fp = B r ■ (area) = B r ■ H ■ 2?rr, where H is 
a semi-thickness of the disk. This flux, within the conducting and differentially rotating Keplerian 
disk, will be wrapped up into a toroidal magnetic field within the disk, Bp- This toroidal field will 
be stronger or a multiple of the initial poloidal field depending upon the number of turns and the 
resistive dissipation of the currents. Initially we consider no dissipation so that an initial flux line 
of poloidal magnetic field Bp will be differentially wrapped n times around the axis leading to an 
enhanced Bt- Let us introduce the number of differential turns, n, that occurs at a radial distance, 
r, during time, t, as 

2vrn = -t ■ r— — , (4) 
dr 

where SIk = (GM/r 3 ) 1//2 is the Keplerian angular velocity of the disk. We consider the toroidal 
flux, Ft, in the azimuthal direction and within the half thickness of the disk, H: Fp = rHBp. 
Then, the increment of this toroidal flux added to the original Ft per dr and per dn differential 
turns becomes 

dF T = dB T ■ Hdr = r ■ B r ■ dt^-^- ■ Hdr = -2ir ■ dn ■ B r ■ Hdr. (5) 

dr 

Since dFT = dBT ■ Hdr, this is equivalent to dBT = —2ndn ■ B r . 

If we integrate over dn and integrate over dr to give the change in toroidal flux, AFt, per 
revolution we obtain an estimate 

AFt ~ -2vrr • n ■ H ■ B r = -n ■ F P . (6) 

The poloidal flux, Fp, in turn is derived from the toroidal flux by the helicity, h, of the plumes driven 
by star-disk collisions. Each plume lifts a loop of toroidal flux with cross section dA = H-Rg^ ~ H 2 , 
Fig. [Q where R s hk is the radius of the shock produced by a star (see paper I for details). The small 
distortion ~ (1 + H/r) from the circular cross section is neglected. When this unit of area or flux 
is rotated ir/2 radians into the poloidal direction with an efficiency of a single plume, a p i ume , it 
creates an equal unit of poloidal flux dFppi ume — — cvpiume ■ Bt ■ H 2 inside the disk. The top parts 
of the field loops created by the plumes are rising quickly because of the low density of plasma in 
the corona: even a very weak magnetic field can overcome the gravity acting on a rarefied gas and 
the expansion will be at a relativistic Alfven speed va % c. Also, the shearing of the top part of the 
loops is small, so the toroidal field produced is also small. The rate of this removal of the toroidal 
and poloidal fields from above the disk to the magnetized jet is higher than the diffusion feeding 
from the inside of the disk, if n / H <C c. We know that tl^H 2 /v ~ 1- Therefore, rj/H ~ and, 
indeed, we have QkH <C c. In this approximation, only the evolution of Fp and Ft inside the disk 
defines the dynamo and it is separate from the evolution of the magnetic fields in the corona of the 



-12- 



disk. Each such unit of flux dFp^ p i ume is only 2H in length and in order to create or to affect a flux 
tube of length r, poloidal, or 2nr, toroidal, requires an aligned sum of increments r/2H in number 
in the poloidal direction and 2irr/2H in number in the toroidal direction. (Each plume also creates 
an increment, a pair of equal and opposite vertical fluxes, dFz — ±dFp tP i ume , which exponentially 
decrease to near zero regardless of overlapping plumes.) 



3.3. Plume Coverage 

The fractional area of one side of the disk inside radius r covered by plumes at any one time, 
q <r , can be estimated as q <r ~ N(< r) • H 2 /r 2 because each of total iV(< r) stars (equation 2 in 
paper I) with impact radii inside a given radius r crosses the disk two times in approximately one 
Keplerian period 2jf(r) = 2tt /f2x( r )- Each such crossing produces one plume of radius ~ H on 
each side of the disk. 

Each plume exists for the time Tk/2 before falling back to the disk surface. Then in the spirit 
of a flux rotation explanation of the dynamo we evaluate the fractional contribution to the total 
poloidal flux by each rotated plume. The plumes occur randomly over the area of the disk, but 
their contribution to the average flux of either poloidal of toroidal is independent of position on the 
surface, because the coherence of the plume rotation ensures an effect proportional to the algebraic 
sum of the number of plumes regardless of their location. Therefore we can rearrange, gedanken- 
wise, the location of the plumes over the disk without affecting the result. We therefore rearrange 
and align a fraction of the plumes, r/2H in number, to create a single, continuous poloidal flux 
tube of poloidal flux, B p i ume H 2 , and length r where B p i ume = a p i urne BT- We have enough plumes, 
iV(< r), to create N tu bes = N(< r)/(r/2H) such flux tubes. These poloidal flux tubes of width 2H 
then collectively cover a fraction in azimuth or a sector of angular width of 

2HN tubes 2H 2 
A(p= — = — ¥ N{<r). (7 

and produce a poloidal flux per half revolution of 

AFp = Ntubes ■ H 2 ■ B p i ume = —Ntubes ■ H 2 ■ OLplume ' &T- (8) 

One can express N tu bes from equation ([7|) as 

TXT 

Ntubes = "^A0. (9) 

We then note that = q <r by construction, and therefore the poloidal flux created by iV(< r) 
plumes per half revolution becomes 

AFp = -Q<r-j^ ■ H 2 ■ a p iumeBT = -nrHq <r a p i ume BT = -q< r ■ ot p i ume -n ■ F T . (10) 

We define a parameter a m as a m = q< r oe p iume to give AFp = — 27ra m • Ft per revolution. Here 
a m becomes an efficiency for the rotation of toroidal flux into poloidal flux per radian of revolution 
and averaged over all the plumes. 
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Then the time derivative considering two plumes per Keplerian period becomes 

dF P 



dt 



Similarly from equation ([6]) we obtain 



-n K ■ a m ■ F T . (11) 



d J^ = -^. 3 -.F P , (12) 
dt 2tt 2 ' v ; 

because one Keplerian revolution, Tk = 2tt/Qk, corresponds ton = 3/2 differential turns according 
to expression (|4]). 

In addition we must consider the fractional flux cancellation of each of these two flux trans- 
formations. This leads to partial cancellation of each orthogonal component by the other and to 
partial self cancellation as well. Finally second order effects, as well as the different dissipation rates 
of the two fluxes must be considered. These effects are usually averaged in mean-field theory, but 
here we consider them separately, because we are concerned with a semi-coherent flow as opposed 
to a turbulent dynamo. 



3.4. Toroidal Multiplication with Losses 

We considered above that the flux is frozen within the disk fluid flow. Now we consider the 
relaxation of this condition by resistive diffusion or reconnection. 

Resistive dissipation of the currents supporting these fields limits the growth of toroidal field. 
Here we consider the saturation of the toroidal multiplication alone with a fixed poloidal field. 
After many turns, this additional toroidal magnetic field reaches a saturation value determined by 
the balance of the multiplication rate with resistive diffusion. The toroidal magnetic field changes 
fastest in the vertical direction on the scale H , therefore the dissipation rate is estimated as 

dFT/dt K .^ Ft _ -jg^-J, (13) 

where magnetic Reynolds number with respect to Keplerian rotation is defined as 

Rm a = (14) 
V 



If we add this loss to the gain of equation (|12p , we then have 



dF T ( 3 F T r 2 , 



dt \ 4tt Rm n H 2 

Thus the toroidal field saturates after KmQH 2 / r 2 turns. In view of equation (|15p the limiting, 
steady state is achieved when 

Ft 3 H 2 
= -— Rmsj— . 

bp 
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This limiting value of the ratio Ft/Fp can be measured separately in the laboratory (in H ~ r geom- 
etry), without the motions producing the a-effect of the complete dynamo, by applying an external, 
fixed poloidal initial field. The same situation may apply to the galaxy if a small residual poloidal 
flux is left over from the initial AGN phase, i.e., this dynamo, an d therefore no further dynamo in 
the galactic disk would be requ i red, even thoug h one likely exists (jRuzmaikin. Sokoloff k, Shukurov 



1988 



Ferriere Schmittl 12000 : iKulsrudl Il999l ) . Furthermore this ratio represents the maximum 



possible toroidal multiplication that should off set losses in the rotation of toroidal flux back into 
poloidal flux for achieving net positive dynamo growth rate. Thus if the toroidal amplification is 
large, the efficiency of rotation of toroidal flux back into poloidal flux, a m , can be small, and the 
dynamo will still be growing. 

To this toroidal multiplication and resistive loss we must add the back reaction effects of the 
helicity or flux rotation mechanism(s), of the "a m -effect". 



3.5. Production of Poloidal Flux and Losses 

In order to calculate a m of equation (fTTj) . we require both q <r and a p i ume . The coverage 
factor, q <r , is more straight forward to estimate from the star disk collision rate and plume size, 
but the efficiency of the helical deformation by a single plume (pQ) is more problematic. The 
simplest and ideal concept of poloidal flux production by a plume is that a plume of radius ~ HJ 2 
rises a distance « 2H above the disk with entrapped flux, dFx jP i U me = Bt • H 2 , rotates this flux 
exactly tt/2 radians, i.e., into the poloidal direction, falls back, merges with the disk matter, and 
releases this now poloidal flux by diffusion or reconnection so that this unit of poloidal flux adds 
in the same direction, i.e. coherently, to Fp. Of course this sequence of rotation, rise, fall, and 
merging of the fluxes will happen episodically and only when averaged, leads to the factor cx m 
such that AFp = —2 Txa m ■ Ft per revolution (equation (|lip ). The associated experimental paper, 



Beckley et al.1 (|2003l ). on laboratory measurements of plume rotation implies that a rising and 
expanding plume, in a rotating frame, indeed rotates a finite angle ~ tt/2 radians before merging 
with the background fluid. In this case the finite angle of rotation occurs for the same reason as 
expected in the accretion disk. In the laboratory case the velocity of the plume relative to the 
velocity of rotation is chosen such that the plume is destroyed or broken up by striking the end 
wall of the apparatus in a chosen finite fraction of a period of rotation. In an accretion disk, as 
pointed out earlier, the plume indeed rises and falls in tt radians with the rotation angle as well as 
the merging with background disk material, both increasing monotonically during this rise and fall 
time. Hence the ideal angle, tt/2 radians, occurs as a result of the product and average of the three 
progressive deformations, but most importantly throughout the entire sequence the incremental 
addition to a m is always positive. Also one should note that with each plume there is an equal 
upward vertical flux, +dFz, as downward vertical flux, —dFz, which presumably averages to zero 
with flux merging. 



As described in equations (jlOl) and (TTTI) the number of plumes adding to the poloidal flux is 
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described by the filling factor, q <r , of the disk by plumes, where q <r = N(< r)H 2 /r 2 . The average 
number of plumes on one side of the disk within a radius r at any given time is approximately equal 
to the number of stars with impact radii inside r, iV(< r), given by expression (2) in paper I. 

H is given by expression (6) in paper I, valid for the inner zone (a) of a standard or Shakura- 
Sunyaev disk. The reader is referred to paper I for th e details, argum ents and caveats of using 
this model of an ac c retion disk originally proposed by IShakural (|1972l ) and further developed in 
Shakura fc Sunvaevi (|1973l ). Here we only remark that we use the expressions for the disk parame- 
ters valid for r < r &, where r a \, is the transition radius between zones (a), (radiation dominated) 
and (b), (particle pressure dominate) generally a few hundreds of r g . As we see from expression (|25p 
below, the dynamo growth rate is maximal at about r a ^. To obtain the number density of stars 
n(r) in the vicinity of CMBH we use analytic and numerical models of the stellar dynamics to 
extrapolate from the observed n at the distance of ~ lpc from CMBH down to few tens of r g . 
Observations typically suggest number densities of the order of 10 5 pc -3 at a distance 1 pc. So we 
write n(lpc) = • 10 5 pc~ 3 . The most notable feature of this distribution of stars is the sharp 



decrease of their density for r less than about lOrt, where rt = 2.1 -10 pc • M, 



1/3 



21r g M 8 



-2/3 



and is the tidal disruption radius for a solar mass star by the tidal forces near the CMBH with 
mass M = Mg • 10 s Mq. This decrease is the effect of physical collisions of stars with each other, 
tidal disruptions by CMBH, and multiple passages of the stars through the accretion disk, which 
grind their orbits into the disk plane and reduce the number of remaining stars not trapped by the 
accretion disk (see paper I for greater details). 

Using the approximate analytical model (2) from paper I we have for r < 10 -2 pc 



<7<r 



1.9 • 1(T 4 • n 5 



for 



l_E_ 

0.1 
r < 10rj 



e 

01 



10r t V ion 
no star-disk collisions), 



-2 



for 10^ < r < 10 pc, 



(16) 



where the factor (1 — ^2>r g /r) coming from the Shakura-Sunyaev model is omitted since r 3> r g . 
Here Ie = L/L Edd is the ratio of the total luminosity of the disk L to the Eddington limit L Edd for 
the CMBH of mass M and e is the fraction of the rest mass energy of the accreted matter, Mc 2 , 
which is radiated away by the disk, L = eMc 2 . The number given by expression (|16p is not large, 
so the probability that any given plume is overlapped with another is small and therefore, on the 
average, each plume will be an individual, isolated event. 

With q <r given by expression (fT6j) the corresponding a m = q< r (Xplume becomes 



10rv 



10r, 



-2 



for 10r( < r < 10 2 pc, 



a m = for r < 10r t (no star-disk collisions). 



(17) 



The system of linear differential equations (jlip and (|12|) has a growing solution 



F T = F Tfi e rt , 



(18) 



-16- 



where 

r = n ir ^ = n^/E^. (19) 

V 47T V 47T 

Similar to the toroidal field, the gradient of the poloidal magnetic field is greatest in the vertical 
direction on the scale H, and therefore the dissipation rate of the poloidal flux is estimated analogous 
to the dissipation rate of the toroidal flux (equation (I13|) ): 



rfF,/* (20) 

where magnetic Reynolds number with respect to the a-deformation is defined as 

„ Qisl 2 , . 

Rm a = (21) 

V 

and / ~ 3H is the height above the disk mid-plane reached by the plume before falling back to the 
disk. In our approximation Rm Q ~ Rmn/ 2 /r 2 , but we keep Rm Q and Rulq separate to evaluate the 
effects of Keplerian and plume motions separately. Adding the resistive dissipation, equation ([20]) . 
to the poloidal gain, equation (fTT|) . results in: 

dFp o j? 12 Qk jp (oo\ 
— = -fl Kam F T - ^—F P . (22) 



The system of linear differential equations (|15p and (|22p has a growing solution of the form (|T8 
but with the growth rate modified as 



2 



'( r 2 P \ 2 3a m V /2 f r 2 I 2 

\H 2 Rm n H 2 Rm a ) vr ) \H 2 Rm n H 2 Rm 



(23) 



We note that in the limit of small resistivity, large magnetic Reynolds numbers, we recover the 
growth rate of equation (|19p . otherwise we note the surprising circumstance that the difference in 
the resistive terms adds to the growth rate whereas, as expected, the sum decreases the growth rate 
as we expect for a purely diffusive resistivity. For positive growth rate, the first term, of course, 
must be greater than the second. In the purely diffusive limit, if one uses that Rm a = Rm^P/r 2 , 
expression ([23|) simplifies to 



f 4i H 2 [\ 4tt H 2 Rm a j V ; 

To the extent that the resistive terms are small and therefore Rm is large and the compensating 
effect of the requirement for merging of newly minted poloidal flux with old poloidal flux is neglected, 

1/2 

then the dynamo growth rate is large, of order Q,k 'dm ■ However, there is no reason to expect that 
the resistivity is purely diffusive, and we expect that tearing mode reconnection drives the merging 



of flu x at some near constant and large Rm ft ~ Rm^ ~ 200 (ISovinec. Finn &: del-Castillo-Negrete 



20011 ) 
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Then from equations (|17p and (|19p the maximum dynamo growth rate occurs between r a b and 
10r t and becomes (for 10 8 M CMBH) 

1/2 



r ~ 7 • io n^t • n/ « p ; ume ' 



0.1 y V0.1/ V 10r t 



-3 

/ r \ 

1 , 

10r t 



for r > 10ri, 



T = for r < 10r t , (25) 
where ri^-j = 2.08 • 10~ 7 s _1 is the Keplerian rotation velocity at lOr^ for a 10 8 M & black hole. The 



exponential rate of multiplication, in view of equation (]25p . maximizes at r = 13.6rt, where 
T max « 1.4 • IO" 9 s- 1 • aj£ me (^) (^) ^ at 13.6r, or 

T ma , « 4 • IO" 2 yr" 1 • «X 2 me (^) _1 at 13.6r t . (26) 

Since the density of stars does not actually drop sharply to at r = lOr^ as in our approximate 
analytical model, the estimate of T max above is approximate and the actual maximum of the growth 
rate is achieved at r somewhat smaller than 13. 6r^. 

One can find the ratio of toroidal to poloidal flux in the growing dynamo mode by substituting 
expression (|23|) for T into equation (fT5"|) : 

Ft = 3/(2tt) 



2 - (27) 

//- fli]) :: ~ .//-Rm„ ) + \j ~~n L + ( H 2 Rm Q ~ J H'2R mct ) 
In the purely diffusive limit, when Rm a = Rm^/ 2 /?" 2 , this ratio simplifies to 

Ft _ B T _ 1 p~ 

F P 27T.BP 2 V vra m ' 1 j 

As in any aw dynamo, the averaged toroidal magnetic flux is much larger than the averaged poloidal 
magnetic flux (recall that a m <1). 

Regardless how small a p i ume , which we believe is ~ 1, the dynamo gain is so large within 
the time of formation of the CMBH, 10 8 years, that saturation will occur early in the history of 
the disk-dynamo regardless how small the initial seed field. The origin of such a seed field, e.g., a 
star, the Biermann battery from decoupling, or primordial fields becomes moot. Nevertheless, for 
completeness, we explore how this gain can be reduced by flux rotated by the plume process such 
that it opposes rather than augments either of the mean toroidal or poloidal fluxes in the above 
estimates of dynamo gain. 



3.6. Flux Compensation by Plumes 



So far we have considered only the positive increment of poloidal flux to the dynamo gain 
by plume rotation of the toroidal flux. This same rotation will rotate the coexisting poloidal flux 
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into opposition with the primary toroidal flux. In addition to the extent that the plumes or any 
other cyclonic motion continues the rotation beyond ir radians, a further reduction in dynamo gain 
occurs because of averaging of this opposed flux. First we consider that Bp/Bt = —(a m /3ir) 1 / 2 <C 1 
according to expression (|28p . 

The deformation leading to the rotation of an increment of toroidal flux into an increment of 
poloidal flux by the rotation of the plumes leads to a similar fraction of poloidal flux being rotated 
such as to oppose the toroidal flux. However, since only a small fraction is rotated by plumes as 
opposed by the shear of rotation, the negative effect on the toroidal flux is small, ~ a m . Similarly 
a rotation by ir radians causes a decrement of both the toroidal and poloidal fluxes to oppose 
themselves, so that the fraction of flux rotated ir radians must be small for high gain. The fraction 
of flux rotated 3n/2 radians must be even smaller for high gain, because a rotation of 3ir/2 radians 
causes the larger toroidal flux to oppose the much smaller poloidal flux even though a small positive 
effect can occur when the same rotation causes the smaller poloidal flux to add to the toroidal. 
This assumes that the plumes are circular in cross section, so that the cross sectional areas for the 
radial and toroidal fluxes are the same. The distortion of the plume cross section by differential 
rotation in ir/2 radians of rotation is similarly small, A(J) ~ (l/A)(H/r). 

We designate these fractions by t±, T2, 73 for the fraction of flux rotated 7r/2, 2tt/2, 3ir/2 radi- 
ans. In general we consider t\ ~ r 2 3> T3, otherwise we do not expect positive gain. One can think 
of these coefficients as correlation coefficients of the decaying plume rotation: a plume undergoes 
little rotation beyond ir radians, when it falls back to the disk. All plumes are considered to be- 
have the same so that these coefficients remain constant and therefore do not describe turbulence. 
These partial flux cancellations all reduce the dynamo gain. The equations (|15p and (|22p are then 
extended to become 

dFx ( 3 Ft r 2 

Mr —r-Fp ~ ^ TFo ~ 0! m ■ n ■ F P - a m ■ t 2 ■ F T + a m ■ t 3 ■ F P 



dt \ 4tt Rm n H 2 

Mr (-am - n- F T - 77-^—775" ~ a m ■ t 2 ■ F P + a m ■ t 3 ■ F T ) . (29) 



dt V Rm a H 2 

Here we have introduced a m • t\ in the first term of the Fp equation where tacitly we had assumed 
T\ = 1 before. Terms with T3 in both equations (|29p are small compared to the terms with T\ 
because T\ 3> T3. Solving system (|29p for exponentially growing solutions we find a generalization 
of expression (|23|) for the growth rate T: 



Ok 
2 



2 , „ v \ 1/2 



/ r 2 I 2 V 

V^^ + ^Rn^ + 2amT2 J. 




(30) 



One can see that the effects of considering finite r 2 and T3 both act to reduce the growth rate of the 
dynamo from the one given by expression (|23|) . Specifically, if T3 > T\ (recall that a m <C 1) then 
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the second term in the sum under the square root in (|30p becomes negative, and T cannot have 
positive real part. This means that the dynamo is impossible for T3 > T\. In fact, T3 enters only 
in combination (n — T3) and effectively reduces the value of t\. This implies that the plume must 
terminate its contribution to flux rotation by tt radians, but this is expected on general grounds 
because by this time the plume matter will have fallen back to and be merged with the disk. 

The effect of finite T2 on the growth rate of the dynamo is much weaker than the effect of 
finite T3. The leading positive contribution to T comes from the second term in the sum under the 
square root and is oc (flmTi) 1 / 2 . The negative contribution of ti is — 2a m T2 term. Since a m <C 1 
we see that this negative contribution will be always small compared to the positive contribution 
for any T2 ~ T\ ~ 1. Thus with T3 assumed small, we expect to recover the very large growth rate, 
r ~ 0.04 per revolution, of expression (|26l) . 



4. Mean Field Theory for the Star Disk Collision Driven Dynamo 



The mean field approach to the problem of generation of t he large scale magnetic fields by the 
motions of the fluid with random component was developed in Steenbeck. Krause &: Radlerl (19661) 
and later was widely used for all possible astrophysical and geophysical applica tions ( Moffatt Il978l ; 



Krause fc Radler 



1980 



Ruzmaikin. Sokoloff Shukurov 



1988 



Kulsrudlll999l ). The basic idea of 



the mean field approach is to average the equations for the evolution of the magnetic field over the 
small scale motions of the conducting liquid. Such small scale motions can be either a collection 
of waves with random phases, or turbulent pulsations, or randomly occurring jets or plumes with 
the sizes considerably smaller than the scale of the whole system. Formal application of the mean 
field theory to the star-disk collision dynamo provides one more mean of justification that such a 
dynamo is operational. 

The number of plumes produced by star-disk collisions is large. At any given moment of 
time there exist ~ 10 4 plumes inside r ~ 10 -2 pc (see Paper I). The radius of each plume is 
r p ~ H ~ 3.7 • 10~ 3 r at r < as shown in section 4.2 of paper I and equation (8) of paper I. 
Therefore, the distance between neighboring plumes is ~ 10~ 2 r and the radial and azimuthal sizes of 
the plumes cannot exceed ~ 10 _2 r without overlapping each other. This condition is well satisfied 
with q <r given by expression (1161) . The magnetic field on the scale of the order of r will be the 
average over many individual plumes. The occurrences of plumes are statistically independent but 
each plume can be considered nearly identical to any other, because the star velocities at any given 
radius are about the same. However, to the extent that the star sizes vary, the energy input to 
each plume will vary accordingly and therefore the size of plumes could be considered as a random 
noise process, but the spectral range is limited. It is attractive to apply mean field theory for the 
generation of the large scale magnetic field by plumes. The averaging over the patches of the disk 
surface exceeding the size of individual plumes is well justified. The averaging over the vertical 
direction is more problematic, since the sign of helicity produced by plumes exactly reverses above 
and below the disk midplane. In addition the typical size of a plume is of the same order as 
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the vertical scale of the change of helicity. Still we explore the results of the application of mean 
field theory equations for the excitation of the global large scale field and attempt to identify the 
departure points of mean field theory from the more coherent flux rotation analysis in this section. 

The random motions induced by the star-disk collisions are clearly statistically anisotropic 
due to the existence of a preferred direction perpendicular to the disk plane as well as a preferred 
direction of rotation on either side of the disk. Still using isotropic expressions for the equations of 
the mean field theory provides so much simplifications that for our purpose of obtaining a proof of 
principle estimate as well as a comparison to flux rotation theory, we will use isotropic equations 
of the mean field theory. The mean electromotive force is given by 



v' x B' = aB - (3V x B, (31) 

where v' are velocities of small scale motions, B' is a small scale field and the bar means averaging 
over small scales (the distances between individual plumes and sizes of the plumes in our case). 
The expressions for the coefficients a and (3 are 

« = "^<v'-(Vx v')>, (32) 

/3 = ^<v' 2 >. (33) 

Here r is the time of the decorrelation of the Lagrangian velocities, i.e. the time of the "memory" 
of a fluid particle about the past history of its velocity. The <> brackets denote averaging over 
the statistical ensemble and in practice can usually be replaced by the averaging over the volume 
larger than the typical scale of the random flow, v', but smaller than the scale of the change of the 
statistical properties of v' and smaller than any large scale of the variability of the mean flow and 
mean magnetic field. 

If the mean large scale flow and large scale magnetic fields are axisymmetric, then one needs 
to solve the following system of equations for the evolution of mean axisy mmetric magnetic field in 



cylindrical coordinates r, </>, z (corresponding unit vectors are e r , e^, e 2 ) ([Roberts fc Sowardl ll992) 
OA 1 / 1 \ 

_ + _vp • V(rA) = ((3 + r] )(v 2 A-^A)+ aB+, (34) 



as, 

dt 



+ rvp • V (l-B^j = rB P • VO + (J3 + rj) (v 2 B^ - ^B^j 



-a I V 2 A - \a) - -Vq • V(rA). (35) 



Here A is related to the poloidal magnetic flux Fp as Fp = 2irrA, B^ is the toroidal magnetic field, 
Bp is the poloidal magnetic field, vp is the poloidal velocity field, and O = Q(r,z) is the angular 
velocity of differential rotation. The quantity A is also a (^-component of a vector potential of the 
mean magnetic field and Bp = V x (Ae^), where is a unit vector in toroidal direction. 

Averaging over the statistical ensemble <> in equations (|32p and (|33|) is replaced by averaging 
over many neighboring plumes. The correlation time r is approximately half of the Keplerian 
period, r = Tk/2. 
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In fact, a and /3 are tensors because of the statistical anisotropy of the plumes. The generali- 
sation of exp ression (|3ip for the mean electromotive force, including effects of anisotropy, is (e.g., 
Moffattlll978l ): v' x B'j = a^B^ — ftijkdBj / 'dxk, where there is a summation over repeated indices. 
In the limit of the qw dynamo, when Rm a <C Rm^, the most important term in the expression 
for the mean electromotive force is a^B^. This term describes the conversion of the toroidal to 
poloidal magnetic flux using the language of the mean field theory. It is analogous to the term on 
the right hand side of equation (fTTI) describing the production of the poloidal flux in the language 
of the flux rotation dynamo. It is this term, which determines the growth rate of the aco dynamo. 
The generation of the mean field in anisotropic random medium is possible for the mean kinetic 
helicity, (v' • (V x y')), e qual to zero, but for non- vanishing compo ne nts of th e tensor fa- (V x v')j.) 
fcrause Radlerl Il98cl : iMolchanov. Ruzmaikin. fc Sokoloflj Il983h . IFerriere! (jl993al lbl. Eiish per- 
formed detailed calculations of the a and (5 tensors resulting from the plume-like motions of gas 
in differentially rotating Galactic disk caused by randomly placed supernovae explosions. These 
motions have some limited similarity to the plumes considered in the present work in that they also 
result in the conversion of the toroidal to the poloidal magnetic flux and are an isotropic due to the 
vertic al density gradient in the Galaxy. Subsequently, these results were used by IFerriere &; Schmitt 
(12000 ) to calculate kinematic anisotropic auj dynamos. In the present work, a and (3 effects are 
assumed to be isotropic. 



We now estimate the magnitude of the coefficients a and (3 in equation (|35p . The half thickness 
of the slab with the helicity produced by plumes is about the vertical extent of a plume, I. We 
assume the dependence of a on z as a = cxqz/1, where ao is a characteristic value of helicity 
which can vary with the radius r. This assumption for a satisfies symmetry requirement that 
a(—z) = —a{z) while exact knowledge of the dependence of a on z is beyond our accuracy. We 
assume that / > H and that a = a^z/l in the whole region —l<z<l, i.e. we neglect the 
fact, that helicity is almost zero inside the disk for —H<z<H. We also assume the turbulent 
magnetic diffusivity, (3, to be uniform over —l<z<l. The fact that the maximum height of the 
plume is I means that the characteristic vertical velocity of the plasma in the plume is v' z ~ vkI/t. 
We assume that the characteristic velocity of the sideways expansion of the plume is v' s ss v' z /2. 
Then, by the time Tk/2 the plume expands to ~ 1/2 in horizontal dimension (we neglect the fact 
that the shape of the plume becomes elliptical). We estimate V x v' ~ — 2f^e z , and therefore 
v' • (V x v') « -2v' z n K = -2in 2 K . Similarly v' 2 = 2vf + v' z 2 ss (3/2)v 2 K (l 2 /r 2 ) for the plume. Let 
us introduce the filling factor q = q(r) equal to the fraction of the surface of the one side of the 
disk covered by plumes. Then averaging, <>, is reduced to the multiplication of the values for one 
plume by q. From expression (|32|) and the above estimate of v' • (V x v') we have 



27r / o 
a = — ■ I ■ Uk ■ q, 

and from expression (|33|) and the above estimate of v' 2 we have 

(3 = l.n K .i 2 . q . 



(36) 



(37) 
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Our estimate of (3 coincides with the estimate of the characteristic value of f3 for an ensemble 
of supernovae explosions occurring at the midplane of the Galaxy considered by iFerrierd ([1993b) 
(formula [351 in that work) . The numerical coefficient in our estimate of (3 is slightly different from 
Ferrierel (jl993bh . 



The dynamo activity is present inside the thin layer with thickness I <C r. This situation is the 
same as for the traditional model of the auo Galactic dynamo. We can use the extensive theory of 
the auj dynamo in thin disks developed in the conn ection with t he Galactic dynamo. An extensive 



treatm ent o f aco Galactic dynamo can be found in lStixl (|1975l ). IZeldovich. Ruzmaikin. &; Sokolofl 



(|1983l ) , and iRuzmaikin. Sokoloff fc Shukurovl (| 19881 ) . One looks for the solution of equations ([31 
and (|35p in the aw limit when Rm a <C Rmj). Since the thickness of the disk, 2H, is small, one 
can neglect radial derivatives of the magnetic field compared to the z-derivatives. In this way 
the problem becomes local with the eigenfrequency of the dynamo determined by solving the one 
dimensional eigenvalue problem in z-direction. This local approximation is similar to the local 
approximation u sed in Appendix A to derive the vertica l structure of the accretion disk. We will 
use results from IRuzmaikin. Sokoloff fc Shukurovl (|1988l ) and replace their parameters with ours. 
The important parameter is the dynamo number 

a l 3 TrQ 2 K ql 4 



D 



dr {j3 + r)f 



(r, + m K Pq) 



(38) 



The D is negative for anticyclonic vortices and d^lx/dr < 0. 



The density of particles in equilibrium non-magnetized disk falls off with z precipitously: oc 
exp(— z 2 /H 2 ) when the gas pressure dominates and even steeper when radiation pressure dominates 
(jShakura &; Sunvaevlll973l ). This means that even a small magnetic field will have a significant 
influence on the dynamics of the disk corona. Thus, the kinematic dynamo approximation does 
not work in the disk corona. There the force-free approximation V x B = AB describes the 
magnetic field evolution at \z\ > I. In the particular case A = the for ce-free magnetic field 



satisfies the vacuum equation V x B = 0. iReves-Ruiz fc Stepinski I (|1999l ) investigated the ato 



turbulent dynamo in accretion disks with linear force-free coronae. They match axisymmetric 
solutions of the dynamo equations (fM|) and (|35|) inside the disk to the solutions with constant 
A of a force-free equation V x B = AB outside the disk. They find that the results for the 
dynamo eigenvalues and dynamo eigenmodes do not change significantly with the value of A. The 
a-quenched saturated mode also depends weekly on A. Thus, in order to obtain estimates for the 
star-disk collisions driven dynamo we can assume that A = and the magnetic fields obey the 
vacuum condition V x B = outside the disk. Note, h owever, that some of the poloidal magnetic 
field lines obtained in IReves-Ruiz &: Stepinski I (|1999l ) have inclination angles to the surface of 
the accretion disk le ss than 60° . This means that MHD outflow should start along these poloidal 
magnetic field lines (jBlandford &: Pavndll982l ). The presence of the MHD outflow would make the 
force-free approximation invalid. However, these field lines, although radial initially, after many 
turns become wrapped up into a force-free helix where the radial magnetic field becomes smaller 
than either the external poloidal or toroidal fields. Both these external fields, in turn are smaller 
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than the toroidal field inside the disk (|Li et al.ll2001al ). Since the magnetic field inside the disk is 
much stronger than outside the disk, the boundary condition at the top of the plume zone, z = ±Z, 
can be approximated as on the boundary with the vacuum: = and B r = 0. 

The eigenvalue pr oblem for the am dynamo in the thin s lab —l(r) < z < l(r) with the vac- 
uum outside the slab (jRuzmaikin. Sokoloff fe Shukurovi (| 19881 )) can be reduced to solving a one- 
dimensional eigenvalue problem in the z-coordinate. In this way, the local growth rate of the 
dynamo T(r) is obtained. The growth rate of the global mode T is very close to the maximum 
value of T m = r(r m ) over the disk radius. The corresponding eigenmode is localized in the ring of 
the disk near radius r m . The characteristic radial width of the eigenmode for the dynamo num- 
bers, that do not much exceed the threshold limit, is ~ (Z^) 1 / 2 ( Ruzmaikin. Sokoloff &: Shukurov 
19881 ). The most easily excited mode of the dynamo has quadrupole symmetry and is steady. The 
excitation condition of this mo st easily excited mode is D < — 7r 4 /16 for the vertical dependence 
of the a-coefficient a = o.qz/1 (|Ruzmaikin. Sokoloff fc Shukurovi Il98a ). The excitation condition 
varies somewhat depending on the choice of the profile of the a-coefficient but is of the same order 
as for the linear profile of a. The growth rate of the most easily excited steady state quadrupole 
mode not far from the excitation threshold is 



P- 




ir 2 



(39) 



The growth rate for large dynamo numbers, \D\ 7r 4 /16, or for small rj is 



0.3 ^ 2 y-n\D\ 



0.3 • Qk^-^/Q- 



(40) 



This differs from Eq. (|T9|) by a negligible factor, ~ 0.35, for a p [ ume = 1, in view of the many 
approximations. We therefore conclude that mean field dynamo theory results in a similar growth 
rate to that predicted by the flux rotation analysis. In either case the growth is so rapid in view 
of Eq. (|26p that nearly the entire history of the accretion disk dynamo will be dominated by the 
near steady state saturated conditions. Unfortunately this steady state is beyond the scope of the 
present paper where instead we feel satisfied in demonstrating an understanding of the dynamo 
gain using a flux rotation model, a mean field theory, and numerical simulations. 

We see that the filling factor q(r) is crucial for the mean field dynamo. Let us estimate q(r). 
The cross section area of the plume is nr 2 ~ ttH 2 , the number of plumes present at any moment 
of time on one side of the disk is 2 • nt>/4 • Tk/2. Therefore, one has 

q = —2—ttH 2 . 
H 4 2 

Using expression (16) of paper I for the flux of stars, nv/4, and expression (A5) of paper I for the 
disk half-thickness, one obtains 

2 



q = 1.52 • 10" 



n 5 



q = for r < 10rj. 



(41) 
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The ratio of the toroidal to the poloidal or radial magnetic field in the growing mode and inside 
the volume occupied by plumes is 



~Bp 



\D\ l l 2 



Using expression (|41|) for the value of q one has 



b~^ 63 ^ VI^ 



-1/2 



nq 



0.1 



-1 



e 

or 



As in all aw dynamos, the generated toroidal field is larger than the poloidal field. However, 
the toroidal field in the vacuum outside the region of dynamo activity vanishes, because the normal 
component of the current at the vacuum boundary must be zero. If there is conductivity, as we ex- 
pect, and therefore force-free magnetic field above the pl ume region, then the toroidal m agnetic field 
generated by the dynamo penetrates into this region fc,eves-Ruiz fc Stepinskil Il99fll ). However, 
due to the quadrupole symmetry of the poloidal magnetic field, the toroidal field in the force-free 
corona has the opposite direction from the toroidal field inside the disk. The axial component 
of the magnetic field, B z , is much smaller than the radial component inside the slab occupied by 
plumes, B z (l/r)B r . However, the radial component of the magnetic field decreases down to the 
value comparable to B z at \z\ = I. The quadrupole poloidal field in the corona is weaker than the 
poloidal magnetic field inside the disk by the factor I jr. The structure of the force- free corona above 
the dynamo generation region cannot be determined without further knowledge about boundary 
conditions at the outer boundaries of the force-free region or physical processes, which limit the 
applicability of force- free ideal MHD approximation in the corona (i.e., fast reconnection of mag- 
ne tic fields). If one requires that the magnetic field in the force- free region vanishes for \z\ ^> I, 



as iReves-Ruiz &: Stepinski I (|1999l ) assume, then, the toroidal magnetic field is comparable to the 
poloidal field in the corona. In this case, the toroidal magnetic field in the force-free corona is much 
smaller than the toroidal magnetic field inside the disk, and so we neglect it in the simulations. In 
the actual case of the black hole accretion disk dynamo, we expect the coronal field to be force-free 
and to progressively remove t h e flux and magnetic energy generated by the dynamo in a force-free 
helix as described in iLi et al.1 (|2001al ) where the field strength, as discussed above, is of the order 
of the poloidal field. 



5. The Dynamo Equations and Numerical Method 

Because of limited numerical resolution and limited computing time we cannot attempt to 
directly simulate the dynamo problem for the real astrophysical parameters. Three dimensional 
simulations of just one star passage through the accretion disk is already quite a challenge for 
computational gas dynamics. Even if we assume that we know the velocity field for a single star- 
disk collision and treat only the kinematic dynamo problem, the existence of ~ 10 4 plumes, the 
necessity of good resolution in the space between and above the plumes, and long evolution times 
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required by the dynamo problem make the direct computations very difficult and demanding of 
major computer resources. Numerical modeling done in this work illustrates and proves essential 
features of the star-disk collisions dynamo described above. We simulate the kinematic dynamo 
with only a few plumes present and adopt a simplified flow model for individual plumes. Then, we 
compare the numerical growth rate and magnetic field structure to the predictions of flux rotation 
and mean field theories extrapolated to a small number of plumes. Qualitative agreement between 
all three approaches in the limit of only a few plumes is observed. 



5.1. Basic Equations 

We have computed order of magnitude estimates of the growth rate and threshold for the 
dynamo by direct numerical simulations. For that purpose we have written a 3D kinematic dynamo 
code evolving the vector potential A of the magnetic field in a given velocity field v and with resistive 
diffusion. The code is written in cylindrical geometry. We start with the equations describing the 
evolution of fields in nonrelativistic quasineutral plasmas. 

V • B = 0, (42) 

if = -V X E, (43) 
Air 

V x B = — j, (44) 

c 

j = a (E + iv x , (45) 

where a is the conductivity of the plasma. Because we are considering the kinematic dynamo, v 
is specified and the momentum equation is ignored. Substituting the expression for the current j 
from the equation (|44p into Ohm's law, equation (|45|) ; and introducing a coefficient of magnetic 

c 2 

diffusivity rj as rj = we obtain Ohm's law in the form 

Ana 

E + -vxB = -VxB. (46) 
c c 

subject to the constraint (l4"2j) . 

The conventional and widely accepted way of writing and solving the kinematic MHD equations 
(MHD without the hydrodynamical part) is to obtain a single equation for the evolution of the 
magnetic field. Substitution of the electric field E from the equation (|46p into Faraday's law, 
equation (|43|) . results in 

— = -V x (r/V x B) + V x (v x B). (47) 
Introducing the vector potential A with 



B = V x A, 
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equation ([37]) takes the form 

OA 

— +r?VxVxA-vx(VxA) + cV</? = 0, (48) 

where tp is the scalar potential; no gauge has been chosen. Any solution of equation (|48p satisfying 
the boundary and initial conditions for the magnetic field should give a physical result for the 
evolution of the magnetic field. The equation (|48|) has the same second order in space derivatives 
as equation ([37]) for the evolution of the magnetic field. 

The gauge freedom can be used to simplify the procedure for solving equation (|48|) . The 
scalar potential ip may be chosen to be an arbitrary function by an appropriate choice of gauge 
transformation. For instance, one can choose to set (p = 0, in which case the remaining equation 
for A takes the form 

dA 

— + r/V x (V x A) - v x (V x A) = 0. (49) 

The boundary conditions for A should be consistent with the gauge chosen. In principle, equa- 
tion ()49p requires three separate boundary conditions for the components of A. This number is 
the same as the number of boundary conditions required to solve the equation for the evolution of 
the magnetic field (j4_7J). Note, however, that there is still a freedom to add Vx to A and there- 
fore to the boundary conditions for A, where x is an arbitrary time independent function, and 
still preserve the gauge condition ip = 0. Although any arbitrary initialization of A satisfying the 
boundary conditions can be allowed, many initializations would result in the same magnetic field 
B. Initializing eq. (I47p V • B = is formally required. We have the following requirements for the 
boundary and initial conditions for A: 

1. There must be boundary and initial conditions on all three components of A. 

2. Boundary and initial conditions should be consistent with the gauge used. 

3. The physical boundary conditions and the initial conditions for the magnetic and electric 
fields (or any other quantities) specific to a particular problem must be satisfied. 

The last requirement means that the physical boundary conditions must be derivable from the 
boundary conditions equations imposed on A. The reverse is not necessarily true, i.e. for one 
specific physical boundary conditions there may be many possible boundary conditions for A. 
The situation with the boundary conditions for A is analogous to the situation with the initial 
conditions for A. With this specification of initial and boundary conditions, the curl of the solution 
to equation ([35]) will be equal to the solution of equation ([37]) . 

If one chooses to evolve the magnetic field directly, then in addition to the equation of evolu- 
tion (]47p the magnetic field must obey the constraint V • B = 0, which should be specified as an 
initial condition. Although it follows from (|47|) that, once initialized to zero, V • B will be kept 
equal to zero, the numerical methods used to solve (|47j) introduce discretization errors, which after 
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a sufficient time can accumulate so that V • B is no longer zero (e.g., lLau Finnlll993l ). Special 
procedures are employed in the codes to deal with this problem such as "divergence cleaning". 
However, in the case of the evolution of the vector potential there are three equations (|48p to solve, 
while there are four dynamic variables in them (i.e. three components of A and one scalar function 
<p). Therefore, one can utilize this one extra degree of freedom in choosing ip for a suitable gauge 
constraint without actually imposing any constraints on three components of A. This will allow us 
to have freedom to choose the gauge and at the same time will not introduce the necessity of taking 
special measures in order to ensure that the gauge will be kept correctly throughout the computa- 
tion. The magnetic field is than obtained by taking curl of A. This way VB vanishes automatically 
within the discretization error associated with approximating the curl by finite differencing. 

In the simulations presented in this work we used the following gauge 

c(^-v-A + t?V-A = (50) 

One can show that for this gauge the basic equation (f4~8|) reduces to 



d A dv k 

— = -A k ^- - (v • V)A + t?V 2 A + (V • A)Vr?. (51) 

We choose the gauge (|50j) because the resulting equation for A has similarity with the equation 
for the advection of a vector quantity. It has the familiar advection term (v • V)A and diffusion 
term 77V 2 A. The term — ^4 fc ^x corresponds to a stretching term (B • V)v in the equation for the 
advection of the magnetic field. Finally, (V • A)V?? term is associated with the nonuniformity of 
electric conductivity. In this work we will consider the case of r/ = constant only and concentrate 
on the effects of the plasma flow producing the dynamo. Thus this term drops out of the equations. 
Note, that the equation (|5ip is valid both for incompressible and compressible flows. 



Finally, we present equations (I5ip written out in cylindrical coordinate system r, cp, z (corre- 
sponding unit vectors are e r , e^, e 2 ) 

dA r ( r dA r 1 6 dA r z dA r 1 . \ / . r 8v r A6 dv^ 
- ' u r -s— + -v^— - + v z — v&Aa - [A r — + A*— — h 



dt \ dr r dip dz r ) \ dr dr 

AZ dv z \ + rq ( l d ( r dAr \ + 1 8,2 A "' + Q2 AT Ar 2 dArP \ + dr l(y A ) ( 52 ) 
dr J \r dr \ dr J r 2 d(j) 2 dz 2 r 2 r 2 d(f> J dr 

dA* _ ( r dA* v^dA^ z dA^ 1 <t >A r 
dt \ dr r d(j) dz r 

A ^ + A *^ + A ^ + W - U'A + (53) 
r dcp r d(p r ocp r r J 

n d f dA<t>\ + 1 d 2 A<i> + d 2 A'*> A* + 2 dA r \ + 1 ^r ] 

\r dr \ dr J r 2 dcj) 2 dz 2 r 2 r 2 dcp J r dcp 

dA z ( r dA z 1 6 dA z z dA z \ ( dv r tA dv* 
v — h -t^^r- + v z -=— - A r — + A*— — h 



dt \ dr r dcp dz J \ dz dz 

A z —\ (-— (r—\ 1 92 A * ° 2AZ \ 9V (V a) (54) 
dz J ^ \r dr \ dr j r 2 deb 2 dz 2 J dz ' 
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1 d 1 dA^ dA z 

where V • A = - — (rA r ) H — — I — - — . The gauge condition (i50l) takes the form 

r or r 6<p oz 

op = v r A r + v t A <t> + v z A * _ rj (--^-(rA r ) + -^r + . (55) 

\r or r oq> oz J 

Also expressions for the magnetic field components in cylindrical coordinates are 

B r = ld^_dM B , = W_dA 1 ^ B* = l^-(rA*)-l^-. (56) 
r dd> dz dz dr ' r dr r d<b 



5.2. Boundary and Initial Conditions 

Although the use of the vector potential eliminates the problem with the divergence clean- 
ing, the boundary conditions in terms of the vector potential may be somewhat more complicated 
and not so obvious from intuitive physical standpoint than the boundary conditions for magnetic 
fields. In this work we used perfectly conducting boundary conditions at all boundaries of the 
cylinder. There is no general agreement on what boundary conditions are most physica l ly ap - 



propriate for a thick accretion disk dynamo simulations. For example, IStepinski &: Levvl ((1988) 
used vacuum boundary conditions outside some given spherical domain for solving the mean field 
dynamo equations in axial symmetry. Khanna & Camenzind (1996a, 1996b) also considered an 
axisymmetric mean field dynamo in the disk and in the corona surrounding the disk on the Kerr 
background gravitational field of a rotating black hole. They used an artificial boundary condition 
that the magnetic field is normal to the rectangular boundary of their computational domain and 
the poloidal component of the current density vanishes near the boundary. However, the main goal 
of these investigations was to demonstrate that certain types of helicity distributions inside the disk 
produce a dynamo. As soon as the boundary of the numerical domain is extended far enough from 
the region of large helicity and large differential rotation, the influence of the boundary conditions 
on the process of the generation of the magnetic fields far inside from the boundary should be 
small. Since both the Keplerian profile of the angular rotational velocity and the frequency of star- 
disk collisions have increasing values toward the central black hole, the approximation of a distant 
boundary can be applicable to the case of our simulations. Therefore we have chosen a perfectly 
conducting rotating cylindrical boundary as a simple boundary condition prescription. We checked 
that the results of our simulations do not strongly depend on the position of the outer boundary. 

The magnetic field near the rotation axis is strongly influenced by the presence of the black 
hole as well as the general relativistic effects associated with the black hole. Magnetic field lines 
in the region close to the rotation axis have their foot-points on the black hole horizon or in the 
region between the black hole and the inner edge of the accretion disk. Therefore, one should 
expect that this region of the magnetosphere will be also strongly influenced by relativistic effects 
of the black hole. The subject of the influence of the central black hole on the magnetic fields 
produced by the dynamo is a part of the so-called "black hole electrodynamics " theory (e.g., see 



-29- 



the chapter "Electrodynamics of Black Holes" in iFrolov fc Novikovl ()1998l )). Since the number 
density of stars should decrease near the black hole due to their capture by the black hole and due 
to tidal disruption, one should not expect the star-disk collision dynamo to operate effectively in 
this region, where strong relativistic effects are important. Therefore, for the purpose of this work 
we replace the region close to the axis of symmetry by imposing an inner cylindrical boundary (also 
perfectly conducting). This may be adequate to the real astrophysical situation in the coronae of 
the accretion disks, since there is highly conducting plasma there. 

We choose as an initial condition a purely poloidal magnetic field with even symmetry with 
respect to the plane of the disk (see Appendix [A] for definitions and properties of odd and even 
magnetic fields). The field is contained within the computational boundaries such that the normal 
component of the magnetic field is zero on all boundaries. 

Let us consider the perfectly conducting rotating boundaries. There is no magnetic flux pene- 
trating the boundaries. This means that the normal component of the magnetic field must always 
remain zero on the boundary. If the velocity of the boundary is Vb, then the tangential component 

of electric field in the rest frame of the moving boundary E H — Vb x B is also zero on the boundary. 

c 

If Vb and B are both tangential at the boundary, then this implies that the tangential component of 
E is also zero there. This then implies that we can chose the <p and the tangential components of A 
to be zero on the boundary. Then from expression (|50p and the vanishing of the normal component 
of v on the boundary, we conclude that we must have V • A = there. Specifically we have 

1 d 

- — (rA r ) = 0, = 0, A z = on the r = constant boundary (57) 
r or 

and 

dA z 

A r = 0, A^ = 0, — — = on the z = constant boundary. (58) 
oz 

This forms a complete set of three boundary conditions for three components of A on each boundary, 
which are compatible both with the physical requirements for fields on a perfectly conducting 
boundary and the gauge condition (j50j) . One can also see that the equations (f57H58j) are consistent 
in the corners of the computational domain, i.e., at the intersections of the planes z = constant 
and cylinders r = constant. 



5.3. The Numerical Scheme 

We use the finite differences predictor-corrector scheme to solve equations (|52H54p in cylindrical 
coordinates. For approximating advection and stretching terms we use central differencing, which 
gives second order accuracy in the coordinates. The diffusion term is approximated by the usual 
7 point stencil. Since the numerical method is explicit, it requires the stability condition to be 

satisfied. Let us denote discretization intervals in coordinates and time as Ar, A<p, Az, and At and 

7]At rjAt r]At v r At v^At v z At 

define the quantities s r = sa = , s z = and C r = — — , La = —r-, C z = — — . 

f r z A(p z Az z Ar rA<p Az 
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Then, the stability conditions that we used in our simulations are 

s r + Sj, + s z <^, {Cr + C^ + C-^) 2 <2{s r + s 4 , + s z ). (59) 

One can show that these conditions follow from the local linear stability analysis of the dynamo 
equations (|52H54p . Before doing each new cycle of predictor-corrector calculations we set up the 
value of the time step At. First, we choose some reasonable value of At dictated by the accuracy 
requirements or how frequent we want to get an output measurements from our simulations. Then, 
we decrease the value of At until the first of the conditions in equation ([59]) is satisfied. After 
that we check the second condition in (|59p and see, if it is satisfied. If not, than we decrease At 
further. One can see, that the second condition in equation (|59p will be always satisfied at some 
value of At since the right hand side depends on At linearly while the left hand side depends on 
At quadratically. The first stability criterion is the usual one for the diffusion equation and means 
that the diffusion per single time step propagates no further than through only a single grid cell. 
The second condition is specific for central differences in the advection term and means that the 
distance the magnetic field is advected during one time step At i s less than the distance through 



which the field diffuses per single time step At (e.g.. I Flet ch eri 1 1 9 9 2| ) . In practice, we ensure stability 



by using a safety coefficient of 0.9 in the inequalities (|59|) . 

When coding the boundary conditions ([57H58]) we used a second order one sided difference 
scheme for approximating derivatives. The resulting expressions have been solved for the unknown 
value of the component of A at the point on the boundary. Boundary conditions have been updated 
after both predictor and corrector steps. In the 4> direction seamless periodic boundary conditions 
have been used, i.e. we make the first and the last grid points in the cj) direction identical and 
corresponding to <j) = and <f) = 2tt and use the same difference scheme as for other values of <fi to 
update these points. Also we used the same seamless treatment of lines = and = 2tt at the 
radial cylindrical boundaries and at the top and bottom boundaries. 

The code is able to treat both the domains with an inner radial boundary and the domains 
including the symmetry axis. In the latter case, there is a singularity of the grid at r = 0, namely, 
all grid points having r = and all values of <j) from to 2ir coincide. One needs a special treatment 
of the grid points at r = ensuring the regularity of Cartesian components A x , A y , A z of A and 
the correct asymptotes for A r , A^ and A z . If the values of the Cartesian components at r — > are 
Ag, Aq, Ag, then the asymptotic behavior of the polar components is A r — > Agcos^ + AqSiikJ), 
A^ — > —Aq sine/; + Aq coscf), A z — > Aq. To impose these asymptotic conditions we first interpolate 
Aq, Aq, and Aq by calculating the average over <p of the Cartesian components of the vector 
potential at grid points situated on a ring with radius Ar. We take this average for Ag, Aq, and 
Aq. Then, we assign the values of the components of A in the cylindrical coordinate system at 
r = according to A r (0) = Agcos^ + A^sin^, A^(0) = -Agsin^ + A^cos^, A z (0) = Ag. This 
finalizes the prescription for the boundary condition at r = 0. When the symmetry axis r = is 
included in the computational region, the code slows considerably because of the small (A^Ar) 
distance between grid points in the (j) direction and, therefore, more restrictive limitations on the 
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time step imposed by the first of the conditions ([59]) . 



5.4. Code Testing 

In the process of writing the code we performed tests for separate parts of the code and, then, 

for the complete code itself. The diffusion part of the code has been tested by reproducing the 

dA 

analytic solution for eigenmodes of the diffusion equation — — = ??V 2 A with A = boundary 

at 

conditions. A variety of different eigennumbers have been tested and decay rates are found to be 
in excellent agreement with analytic expressions. The code preserves the shape of eigenmodes with 
very high accuracy even for a very moderate number of nodes. Coupling between equation (|52[) for 
A r and equation (|53|) for has been tested by evolving nonaxisymmetric eigenmodes. 

The advection part of the code has been tested by computing the advection by the uniform 
flow of the magnetic field of the type B = Bn, where n is a fixed vector of unit length (we made 
a few runs with different directions of n), and the magnitude of the magnetic field B has the 
constant gradient vector V-B = constant perpendicular to n. The current density corresponding to 
such a magnetic field is uniform, and therefore, the magnetic field does not diffuse. The boundary 
condition for this test was set to time-dependent explicit values computed from the known purely 
advective behavior of the field. We observed good agreement with the picture of the pure advection 
of flow. 

We also compared the results for dynamo simulations with the two dimensional flow given 
by our code to the simulations produced by two other 2D kinematic dynamo codes, one evolving 
vector potential and another evolving magnetic field. The latter 2D code has a divergence cleaning 
procedure for VB. The flow was an axisymmetric Beltrami flow with V x v = Av. For the interior 
of the domain < r < R Q and < z < L one can obtain the following analytic solution for the 
Beltrami flow: 

v r = J\ (jn 
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where Jq(x) and JAx) are the Bessel functions, jn is the first root of J\(x) = 0, A| = — \ + — r. 
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The solution can also be written in terms of the flux function ^(r, z): 
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The 3D kinematic code picks up the fastest growing mode of the dynamo. In the case of axisymmet- 
ric flows the nonaxisymmetric modes of the field (oc g* 7 "^ - * 1 ^') with different azimuthal wavenumber 
m evolves separately. The fastest growing mode in our simulations was with m = 1. The growth 
rate and the structure of the m = 1 modes obtained with 3D and 2D codes agrees remarkably 
well. We also studied the convergence with respect to the grid resolution and found that for a 
magnetic Reynolds number R m (defined as the product of maximum velocity and minimum of L 
and R ) of about 200 the simulations converge for the grid resolution of about 41x61x41 in r,(p, 
and z directions respectively. 



6. Results of Numerical Simulations 

6.1. Model of the Flow Field 

We now approximate the flow for our kinematic code from the analysis of the simplified model 
of plumes in Section [3j When describing the results of our numerical simulations, we will use 
dimensionless units with the unit of length equal to the radius at which the star-disk collisions 
occur, and the unit of velocity equal to the Keplerian velocity at that radius. Then, one turn of 
the disk at unit radius takes 2ir dimensionless units of time. The disk is assumed to have constant 
thickness. Its top boundary is at z = zt op and bottom boundary is at z = Zbot- We usually put the 
disk at z = 0, in the middle of computational cylindrical domain, and then, Zbot = —Ztop- However, 
we will preserve separate notations for top and bottom boundaries. For simplicity we assume that 
all star-disk collisions happen at unit radius, but are randomly distributed in azimuthal angle along 
r = 1. Also, a remarkable feature of star-disk collisions is that the numbers of stars crossing the 
disk in both directions are equal on average. We consider two models for the position of star-disk 
collisions addressing this property. In the first model we assume that collisions happen in pairs: at 
each time there are two collisions at r = 1, one with the star going up through the disk and the 
other at the opposite point on the circle r = 1 with the star going down through the disk. Thus, 
at any moment of time the flow is symmetric with respect to the inversion relative to the central 
point of the disk. The second model considers random directions of plumes as well as random 
distribution of plumes over the circle r = 1. In the next section we describe the results obtained 
with both models. 

The plume flow is superimposed onto a background of Keplerian differential rotation occupying 
the whole computational domain = A star-disk collision is simulated by a vertically 

progressing cylinder of radius r p in the corotation frame. The cylinder starts at the bottom of the 
disk located at z = Zbot, penetrates the disk, and rises to a height of h above the disk. At the 
same time the cylinder rotates about its axis opposite to the local Keplerian frame such that the 
cylinder does not rotate about its axis if viewed in laboratory frame (an inertial frame where the 
central black hole is at rest), but the axis corotates with the local Keplerian frame. By the time 
the plume reaches its highest point, n/2 radians of Keplerian rotation, the axis corotates with the 
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Keplerian flow also by tt/2 radians on average. Since the cylinder does not rotate about its axis, 
the relative rotation between the cylinder and Keplerian flow corresponds to an untwisting of tt/2 
radians, when the local frame rotates tt/2 radians as measured at the radius of the axis of the jet. 
The length of the cylinder is progressive with time and its velocity, v pz ~ vk- The vertical velocity 
of the gas inside the cylinder is constant and is equal to v pz . After the time the plume rotates 
by tt/2 it is stopped and the velocity field is restored to be pure Keplerian differential rotation 
everywhere. This very simplified flow field captures the basic features of actual complicated flow 
produced by randomly distributed star-disk collisions. We also feel that elaborating on some of 
the details of the flow field like taking a more realistic distribution of star-disk collision points in r, 
and introducing a weak and distributed downflow, is not warranted at the present initial stage of 
simulations in view of the fact that we do not know other important features of the flow (no actual 
hydrodynamic calculations have yet been performed ). Our model flow and simplified assumptions 
about star-disk collisions, frequency, and distribution capture qualitative features important for the 
excitation and symmetry properties of the dynamo. We feel that all elaborations mentioned above 
as well as accurate simulations of star-disk collision hydrodynamics would not qualitatively change 
our conclusion about the possibility of such a dynamo. 

Since equations (f52H5"i| require spatial derivatives of the velocities, we apply smoothing of 
discontinuities in the flow field described above. Also we introduce smooth switching on and off of 
the plumes in time. For all three components of velocity v k we use the same interpolation rule for 
two plumes 

v k = v k mlSl + v f n2 s 2 + (1 - si - s 2 )v k out . (60) 

Here si(r, (ft, z, t) and S2(r, (ft, z, t) are smoothing functions for plume 1 and 2 correspondingly. Each 
function s is close to 1 in the region of space and time occupied by the plume and is close to in 
the rest of space and during times when the plume is off. Transition from 1 to happens in the 
narrow layer at the boundary of the plume and during the interval of time short compared to the 
characteristic time of the plume rise. v k nl and v\ n2 are velocities of the flow of plumes 1 and 2, v k ut 
is the velocity of the flow outside the regions occupied by the plumes. For spatial derivatives of the 
velocity components, one has from (|60p 

dv k _ ds x ( k k \ ds 2 / k k \ 8v k inl dv k n2 ( ,dv k out 

It is easy to generalize this approach for an arbitrary number of plumes. 

Let us assume that the position of the axis of a cylindrical jet launched upward (in the positive 
direction of the z axis) is at r = r$ and (ft = 4>$. We keep ro = 1 for all plumes and the initial (fto is 
randomly taken between and 2ir. Let us denote this plume as number 1 and the symmetric plume 
going down from the equatorial plane as number 2. Then, after time (t — t p ) from the starting 
moment of the plume t = t p , its position is 



<ft\ = (fto + (t - t p )r £l K o, 



(62) 
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where £Iko = ^E"(^o) is the Keplerian angular rotational velocity at r = ro and in the simulations 
presented in this work, CIko = 1- The position of the axis of the symmetric plume is 



4>2 = <Pl + 



(63) 



The radii of both plumes are r p . The bottom surface of the plume 1 is at z = Zbot, the top surface 
of the plume 1 is at z\ = z^t + v pz (t — t p ), the top surface of the plume 2 is at z = z top , the bottom 



surface of the plume 2 is at Z2 
inside the upward jet is 



Ztop 



"pz 



(t — t p ). Due to symmetry, Z2 = —z\. The velocity field 



v\ = r tt K0 sin(> - <f>i), 
vf = r tl K0 cos((p - 4>i), 



v\ =- r 



pz- 



The velocity field inside the downward jet is 



t>2 = Tq^IkO COS(^) 

v z 2 = -v„ 
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We choose the following interpolation functions 
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Here r[ 2 = r^+r 2 — 2r$r cos(0 — (pi) is the distance from the axis of the plume 1, r 2 2 = r\ + r 2 — 
2rorcos(</> — <p2 ) is the distance from the axis of the plume 2, A is the thickness of the transition 
layer of the functions si and s 2 from their value 1 inside the plume to outside the plume, A <C r p . 
Square root expressions in the z-parts of si and s 2 ensure that the thickness of the transition layer 
in the z direction is never less than A, even just after the plumes are started, when the differences 
(zi — Zbot) and (ztop — z 2 ) are zero. We choose A = 0.01. 



The function S(t) ensures smooth "turning on" and "turning off" of the plumes at prescribed 
moments of time. If the plumes are to be started at t = t p and to be turned off at t = (t^ > t p ), 
then we adopt the following form of the function S(t) 
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where 5t is the length of the transition period. S = corresponds to the flow without plumes, 
and S = 1 corresponds to the flow with plumes. One needs to ensure that 5t < td — t p . We took 
St = (td — t p )/5. The cycles with the cylindrical jets present are interchanged periodically with the 
cycles with the pure Keplerian rotation only. The time between two consequent launchings of the 
plumes is At p and we always have At p > t^ — t p , such that at any time only one pair of plumes is 
present. This eliminates the occurrences of overlapping jets. Note, that during the time td — t p the 
disk makes only about a quarter of the turn. 

Our second model of random directions of the plumes introduces obvious changes into the 
expressions above. Namely, we set S2 = in equations (|60p and (|6ip . and we intermittently use 
either expressions ([63H66]) for the velocity, when the jet is directed upward, or expressions ([67H69]) . 
when the jet is directed downward. We also use the same "switch" function S(t) for both models. 

Finally, let us list the parameters, which are important for the growth of the magnetic field 

r 2 Q K (r) 

in our model: the magnetic diffusivity rj, (or magnetic Reynolds number Rm^ = ), the 

V 

radius of the plumes r p , the frequency of star-disk collisions, At p , the vertical velocity of the plume 
v pz , and the duration of the plumes td — t p . 



6.2. Analytic Solution in the Asymptotic Region 

Because the equations for the evolution of the magnetic field B and vector potential A are 
parabolic, the boundary conditions will always influence the solutions inside the computational 
region. However, the distribution of the frequency of star-disk collisions is concentrated towards 
the center meaning that most of dynamo activity happens in a limited region of space (around 
r = 1 in our dimensionless units). If one is willing to disregard a relatively small a-effect for r>l, 
then the solutions of the field equations in the region r>l can be obtained analytically because 
the flow is just the differential rotation with the Keplerian angular velocity. 

The equations for the evolution of the axisymmetric magnetic field in the presence of only the 
differential rotation are analogous to equations (|34p and (|35p for the evolution of the axisymmetric 
mean field. We can obtain the necessary equations when replacing the mean field by the actual 
field and using the same functions A and for the poloidal magnetic flux and toroidal magnetic 
field. If one sets a = 0, f3 = 0, vp = 0, and £1 = in equations (f3~ij) and ([35]) . the resulting 

equations for axisymmetric magnetic field in a purely rotating flow are 



Equation (|72p is a diffusion equation for the poloidal magnetic field without sources. Its solutions 
are determined by boundary conditions imposed on the poloidal magnetic field. Equation (|73p is 
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a diffusion equation for the toroidal magnetic field but with the source term due to the fi-effect. 
We see that the evolution of poloidal magnetic field is decoupled from the evolution of the toroidal 
magnetic field (unless boundary conditions mix them together). After one knows the solution for 
the poloidal field, one can solve equation (|T3|) to find the toroidal magnetic field. If one looks for 
stationary solutions of equations (|72|) and (|T3|) then the outer boundary condition is very important 
to determine the solution. However, in the case of a dynamo the magnetic field in the dynamo 
domain r ~ 1 grows exponentially. This growing field diffuses into the surrounding conducting 
medium according to equations (j72[) and (173)) . The phenomenon is analogous to the skin layer in 
plasma. The growing magnetic field decreases exponentially outward from the generation region. 
Therefore, if the growth rate is sufficiently high such that the skin depth is smaller than the distance 
to the ideally conducting boundary, the boundary conditions at the boundary do not influence the 
dynamo process. 

We computed an analytic solution of equation (|72l) in the region r > 1 when the magnetic field 
grows exponentially. This solution is presented in Appendix [Bj We have checked with numerical 
simulations of the dynamo that the magnetic field in the zone outside of dynamo activity but inside 
the outer radius of our computational domain is very closely approximated by expressions (|B6j) 
resulting from our analytic solution. We also varied the outer sizes of the outer ideally conducting 
boundaries in our 3-dimensional simulations to verify that the growth rate and the structure of 
the growing magnetic field are insensitive to the placement of the boundaries. It is necessary to 
stress that the simulations are insensitive to the boundary conditions only when the magnetic field 
is exponentially growing: the simulations in the cases of decaying or steady fields do depend on 
how far the ideally conducting boundaries are placed. 



6.3. Simulations of the Dynamo Growth 

The simulation is shown in a sequence of stages. We use dimensionless units described in 
section 16.11 Our computational domain is the space between two cylinders with the inner radius 
R\ = 0.2 and the outer radius R2 = 4, filled with a media having uniform magnetic diffusivity r\. 
The computational space is limited from below by the surface z = —4 and from above by the surface 
z = 4. The total length of the cylindrical volume comprised between surfaces z = —4 and z = 4 is 
8. All boundaries of the computational volume are ideally conducting. There is no magnetic field 
penetrating the boundaries, and the boundary conditions (|57|) and ([58]) are applied. 

An initial quadrupole like field establishes a primarily radial field within the midplane of the 
cylindrical volume, \z\ < 1/3. The initial field is purely poloidal, concentrated toward the inner 
parts of the disk, and is shown in Fig. [2] by arrows. The accretion disk is indicated at Zt Q p — 

1/3, 

Zbot = -1/3. 

Keplerian differential rotation is initiated and generates toroidal field. At the same time 
poloidal field diffuses toward the outer boundary and becomes distributed over the volume more 
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uniformly. The magnetic diffusivity is rj = 0.01, and the magnetic Reynolds number for rotation 
r 2 Q K (r) 

at r = 1 is Rhiq = = 100. With no source term present in equation (|72|) . the poloidal 

rj 

magnetic field will decay away in a purely toroidal flow. The toroidal magnetic field will 
first grow because of the source term r — r—B r in equation (I73h . then reach a saturation value 

~ i?pRmoy(2-7r) determined by the balance between the source term r — r^-B r and the diffusion 

dr 

term rj(V 2 B ( j ) — r~ 2 B$) in equation (l73l) . and finally decay as the poloidal magnetic field B r decays 
and so the source term for also decays (Cowling's theorem). 

Fig. [3] illustrates the poloidal magnetic field obtained after several revolutions at r = 1, Fig. [J] 
shows the contours of toroidal field at the same moment of time as on Fig. [3] The time evolution of 
the fluxes of magnetic field is shown in Fig. We also show the process of winding up the dipole 
like (odd) field in Figs. [6] and [7] (poloidal and toroidal fields) and Fig. [8] (the evolution of fluxes). 
Note, that the toroidal field produced from the initial quadrupole field (and any even symmetry 
field) has the same sign throughout the disk thickness as well as in the space above and below the 
disk. In contrast the toroidal field produced from an initial dipole field (and any odd symmetry 
field) is zero at the equatorial plane and has opposite signs in the upper and lower halfs of the disk 
thickness. 

We now examine how the simulated star-disk collisions (approximated by a flow model de- 
scribed in section RTTj) deform the wound up, toroidal magnetic field and create poloidal field from 
the toroidal field. Figs. [9] and [10] illustrate the action of the rising plume on the poloidal magnetic 
field in a fluid which is at rest. The initial magnetic field here is a quadrupole like field shown in 
Fig. [2] The radius of the inner cylinder is 0.2 and the radius of the plume is 0.2. The velocity 
of the plume is equal to the Keplerian velocity at r = 1 and the plume moves tt/4 radians in the 
(^-direction before it disappears. Fig. [9] is a side view on the plume. Fig. [10] is a view on the plume 
from the top. One can clearly see the lifting of the field lines of the quadrupole field from the 
midplane of the disk by the plume flow. Because the plume flow is strongly compressible near the 
head of the plume it forms a narrow layer of enhanced magnetic field near the top boundary of 
the plume. Magnetic field diffuses inside the plume from this layer. On the top view one can see 
the twisting of magnetic field lines by the unwinding of the flow in the plume. It creates toroidal 
field from the poloidal field. More importantly, Figs. [TT1 [T2l and [13] illustrate the action of the 
same plume on the primarily toroidal magnetic field wound up from the initial quadrupole field 
(as in Fig. Hj). The plume rises through the differentially rotating fluid with the Keplerian profile 
of angular velocity. Fig. [11] is a side view from r-direction, Fig. [12] is a top view from z-direction, 
and Fig. [13] is a side view from ^-direction. Shown by arrows is the flow velocity in the reference 
frame corotating with the base of the plume with the angular velocity at the point of the location 
of the plume, i.e. the value v' = v — Qxo re </>- As with Fig. [9] the side view from the r-direction 
on Fig. [11] shows the lifting up of the toroidal field by the rising plume. One can see from the 
projection viewed from (^-direction that the magnetic field is entrained into the forming a loop of 
poloidal field. The top view clearly shows the twisting of toroidal magnetic field and the creation 
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of poloidal field from the toroidal field, i.e. the a effect. The resulting loop of flux translated and 
rotated from the toroidal plane is shown at the time of maximum jet extension. After that time 
the jet velocities are smoothly set to zero. 

By close examination of the positions of field lines in Figs. [TT j 1121 and [13] one can discover the 
presence of another, more subtle effect: as the bundle of magnetic field lines is rotated and bent by 
the plume, magnetic field lines twist around each other in this bundle. The direction of this twist 
can be observed to be opposite to the direction of the helical twisting associated with the lifting 
and bending of the bundle as a whole. The bundle of magnetic field lines behaves like a ribbon 
when it is bended and curved. The reason for the additio nal opposite twist of the magnet ic field 



lines in this ribbon is the conservation of magnetic helicity (jBlackman &: Brandenburgj|2003l ). This 
small scale twist does not influence our flux rotation and mean field estimates of the kinematic 
stage of the dynamo. 

The problem is continued with the jets or plumes repeated. The model of the flow described in 
section [pTTI is applied. Below we present the results for a representative case for the model with the 
plumes randomly distributed along the circle r = 1 and launched in periodic intervals in random 
directions up and down through the disk. The parameters are the following (in dimensionless units 
introduced in section [6T}) : R\ = 0.2, R 2 = 4, rj = 0.01, r p = 0.3, At p = vr/2 + 0.4, t d -t p = vr/2, 
v pz = 1, Zb t = —1/3, ztop = 1/3 and the centers of plumes are located on the circle r = 1. The 
run is started with the initial field being purely poloidal. The initial poloidal field is the linear 
superposition of odd and even magnetic fields shown in Fig. [6] and Fig. [3] respectively. The exact 
meaning of odd (dipole like) and even (quadrupole like) parity fields is described in Appendix El 
Here we only note, that the total energy of the magnetic field is equal to the sum of energies of odd 
and even components. Odd component contributes 5% of the total energy of the initial field. The 
remaining 95% of the total energy is the energy of the even field. The first plume is launched at 
the moment t = 0.2 after the beginning of the simulation, and the subsequent plumes are launched 
in periodic moments of time with the period At p . This rate of plume launches corresponds to an 
average 2ir/At p = 3.2 plume launches per revolution at r = 1. The simulation is continued until 
time t = 640. By that time the magnetic field grows by ~ 10 orders of magnitude. The resolution 
of our typical dynamo simulation is 41x81x41 nodes in radial, azimuthal and vertical directions 
respectively. Although this resolution seems to be quite modest to resolve the plumes (there are 
typically only about 6x6 nodes to resolve the cross section of a plume) we checked the convergence of 
our simulations by performing trial runs with 61x121x61 resolution. The growth rate of the dynamo 
and the structure of the growing magnetic fields do not change with the increased resolution. We 
also performed trial runs with the larger size of the computational domain: —6 < z < 6 and 
0.2 < r < 6 with 61x121x61 resolution. We did not observe significant changes of the growth rates 
and magnetic field structure of the dynamo when increasing the size of computational domain. The 
reasons for insensitivity to the boundary conditions are described above in section I6.2L 



The time evolution of the total energy of the magnetic field integrated over the computational 
volume is presented in Fig. [J4] as well as the time evolution of the fractions of total energy of 
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odd and even components of the magnetic field. An arbitrary value of the initial magnitude of the 
magnetic field is used. The initial rapid growth of the energy is due to rapid build up of the toroidal 
magnetic field. After a couple of revolutions at r = 1 the dynamo effect overcomes the linear growth 
of the toroidal magnetic field and the growth of the magnetic energy becomes exponential. The 
magnetic field experiences oscillations with the period equal to At p due to the repeated actions 
of single plumes. More significant oscillations of odd and even components of the field occur on 
the time scale of the diffusion over the region of dynamo activity ~ 100. Despite the significant 
variation of the fraction of the odd field, which can become up to 30%, even (quadrupole) field 
dominates. Since the flow does not have symmetry with respect to reflections z — > — z, the odd and 
even components of magnetic field are coupled to each other and grow with the same exponential 
rate. 

The time evolution of fluxes of three components of the magnetic field is shown in Fig. [15] We 
calculate the fluxes of magnetic field through the following three surfaces: the flux of B r through 
the part of cylindrical surface r = 1/2 limited by lines z = 0, z = 4, (j) = 0, and <j) = it/2; the flux of 
through the rectangle in the plane <fi = limited by lines z = 0, z = 4, r = R\, and r = R2; the 
flux of B z through the half of the ring in the plane z = — 2 limited by lines r = R\, r = R2, <fi = 0, 
and <j> = it. Then, we divide each of the three fluxes by the areas of the corresponding surfaces. 
In this way, the values of the magnetic field averaged over the surfaces, < B r >, < B^ >, and 
< B z >, are obtained. The time evolution of the logarithms of absolute values of these averaged 
values of the magnetic field is presented in Fig. [15] All three fluxes grow exponentially (if averaged 
over fluctuations) with the same growth rate V = 0.026. The growth rate of the mean square of the 
magnetic field plotted in Fig. [14] is equal to 2r which is consistent with the growth rate of fluxes. 
The value of < B^ > is larger than the values of poloidal fluxes meaning that the toroidal field is 
predominant in the dynamo, which is also in the agreement with the conclusion from the mean field 
theory. While radial and toroidal fluxes grow monotonically, the flux of the axial magnetic field 
experiences oscillations with exponentially growing amplitude. The z-flux remains zero on average. 
This is due to the fact that both dipole and quadrupole growing magnetic fields have zero z-flux 
through the surface described above. However, the z-flux experiences oscillations due to individual 
plumes creating nonaxisymmetric magnetic field. 

The behavior of dynamo magnetic fields immediately outside of the generation region is espe- 
cially interesting in connection to the magnetic fields in the jets (magnetic helices) and observed 
magnetic field in galactic disks. In Fig. [16] we plotted the fraction of energy of the magnetic field, 
which resides outside of the region of dynamo activity. In particular, we divided the whole com- 
putational domain into two: the inner domain is the region —2 < z < 2 and r < 2, the outer 
domain is the rest of the computational domain with \z\ > 2 or r > 2. Initially, the fraction of the 
outer energy grows because of the diffusion of the initial magnetic field outside the central region 
(compare the poloidal field on Fig. [2] and on Fig. [3]). However, after the dynamo action sets in, the 
skin effect described in section [531 and Appendix [B] occurs. The skin depth of the steady growing 
magnetic field given by equation (|B8|) for rj = 0.01 and V = 0.026 is / s = 0.6. Thus, the outer 
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domain is in the zone of pure diffusion of the magnetic field, where the variations due to individ- 
ual plumes are smoothed out. The average value of the outer fraction of the magnetic energy is 
~ 0.06 of the total magnetic energy. This is roughly consistent with the estimate one can obtain 
from the skin depth analysis of Appendix [B] ~ (0.6/e) 2 ~ 0.05. The field in the outer region is 
predominantly even as well as in the inner region. The time dependence of the fraction of even 
field in Fig. 1 161 follows closely the time dependence of the fraction of the even field in Fig. [T4j Note, 
however, that the curves in Fig. [16] are more smooth than in Fig. [T4j Rapid oscillations of the field 
caused by individual plumes are smoothed out in the diffusion process of the magnetic field into 
the outer region as the exponential decay scale l s becomes shorter for higher oscillatory frequencies 
to' (Appendix [Bj. Only slow variations with the time scale about or longer than the diffusive time 
scale remain present in the outer domain. 

Another diagnostic of our simulation is to calculate the time behavior of the magnetic fluxes 
through the surfaces in the outer part of computational domain. By looking at the time evolution of 
these fluxes we can learn about the time evolution of the magnetic field in the asymptotic diffusion 
region. We calculate magnetic fluxes of radial magnetic field, or equivalently, < B r > through the 
following cylindrical surfaces: radial flux 1 through the part of the surface r = 2 limited by lines 
(ft = 0, (ft = 7r /2, z = —1/3, and z = 1/3; radial flux 2 through the part of the surface r = 3 limited 
by lines (ft = 0, (ft = tt/2, z = —1/3, and z = 1/3; radial flux 3 through the part of the surface r = 3 
limited by lines (ft = 0, (ft = tt/2, z = 2, and z = 4; radial flux 4 through the part of the surface 
r = 3 limited by lines (ft = 0, (ft = tt/2, z = —4, and z = —2. The first two radial fluxes describe the 
evolution of the magnetic field close to the equatorial plane. The third and fourth fluxes describe 
the evolution of the magnetic field in the outer corners of the computational domain. We plot these 
four radial fluxes in Fig. [T7l We calculate three fluxes of the toroidal magnetic field, or equivalently, 
< -E?0 > through the following rectangular areas of the plane (ft = 0: toroidal flux 1 through the 
rectangle limited by lines r = 2, r = 4, z = — 1/3, and z = 1/3; toroidal flux 2 through the rectangle 
limited by lines r = 2, r = 4, z = 3, and z = 4; toroidal flux 3 through the rectangle limited by 
lines r = 2, r = 4, z = —4, and z = —3. We plot these three toroidal fluxes in Fig. [18] We calculate 
two fluxes of the axial magnetic field, or equivalently, < B z > through the following ring-shaped 
surfaces: axial flux 1 through the quarter of the ring in the plane z = 2 limited by the lines (ft = 0, 
(ft = tt/2, r = 3, and r = 4; axial flux 2 through the quarter of the ring in the plane z = —2 limited 
by the lines (ft = 0, (ft = tt/2, r = 3, and r = 4. We plot these two axial fluxes in Fig. [T9l One can 
see that all radial, toroidal and axial fluxes do not change sign during the exponential growth of 
the dynamo (after the time t ~ 100). Therefore, the star-disk collisions dynamo produces steadily 
growing non-oscillating magnetic fields. The signs of the fluxes (not shown in Figs. [TTMl9[) are 
consistent with the quadrupole geometry of the magnetic field in the outer region of the dynamo. 

In Fig. [20] we plotted two vector plots of the poloidal magnetic field at the plane (ft = at the 
final moment of the simulation t = 640: on the top plot the length of arrows is proportional to the 
magnitude of the poloidal magnetic field, on the bottom plot all arrows have unit length and the 
direction of the arrows indicate the direction of the same magnetic field as on the top plot. The 
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concentration of the magnetic field toward the central region with the plumes is clearly visible on 
the top plot. The imaging with arrows picks up only the region of the strong field while the arrows 
outside this region are so short that they cannot be pictured at all. The bottom plot illustrates the 
structure of the poloidal field in the asymptotic outer region. This structure can be described as 
a "shifted quadrupole" implying the presence of a significant dipole component. The toroidal field 
is ~ 20 times stronger than the poloidal. The direction of the toroidal field agrees well with the 
direction of the field produced by the stretching of the poloidal field by the Keplerian differential 
rotation. The structure of the field at different <ft positions is similar to that at <p = 0. The 
nonaxisymmetric variations of the field are most significant at the location of the plumes at r ~ 1 
and quickly decay outwards. Each individual plume perturbs the magnetic field significantly. This 
is also reflected in the oscillations of fluxes in Fig. [15j The three dimensional plot of the dynamo 
magnetic field is presented in Fig. [2T1 Here we plotted only the poloidal component of the magnetic 
field at the two meridional slices, (p = ir/2 and (ft = 3ir/2, in the computational domain. In order 
to smooth out the strong contrast between magnitudes of the magnetic field in the inner and outer 
regions of the computational domain, we plotted a vector field Bp/|i?p| 2//3 . The dominance of the 
quadrupole magnetic field in the outer asymptotic region is obvious from Fig. [2X1 In the central 
region for r ~ 1, the field is strongly perturbed by individual plumes, and the nonaxisymmetric 
field caused by the action of each single plume is visible. Toroidal magnetic field is also strongest 
in the central part of the computational domain. 

Finally, let us compare the predictions of the flux rotation and the mean field theories with 
the results of our numerical simulation. All three predict that the growing magnetic field will be 

quadrupole. The simulation formally corresponds to q « q <r = — - — - = 0.036, H = —z\, t = 

r A 2At p 

1/3, I = Zbot + v pz (t(i — t p ) = 1.24 in dimensionless units of simulation. Using these parameters 
and OLpi urae = 1 in the expression for the growth rate in the flux rotation theory, equation (124p . 
one obtains V = 0.084. For the mean field theory the expression (|3T[) gives = 0.09, the dynamo 
number (equation (f38jl ) is D = —28, and both expressions (f39~|) and (|4"0j) give T « 0.18. This is to 
be compared to numerical growth rate T = 0.026. Both the flux rotation and especially the mean 
field theory growth rates are higher, but all three are within one order of magnitude from each 
other. Such a result is satisfactory because of the far reaching extrapolations of the applicability 
of both flux rotation and mean field theories. 



7. Conclusions 

We believe that by theory and calculation we have demonstrated that a robust aui dynamo is 
likely to occur in conducting accretion disks with a robust source of helicity. The growth rates as 
large as V ~ 0.1 to 0.01 VLk are expected. We have discussed in depth one such source of helicity 
in the accretion disk forming the central massive black holes of most galaxies. This is the almost 
inevitable star-disk collisions that should occur in the dense stellar populations at the center of the 



-42 - 



galaxy. We estimate that this source of helicity is far larger than necessary for the dynamo fields to 
reach saturation in less than the formation time of the black hole. Star-disk collisions should also 
be the most robust source of helicity because the resulting plumes are driven several scale heights 
above the surface of the disk as compared to turbulence were the vertical motions are limited to 
a fraction of a scale height. The advantage of the au: dynamo is that because it produces a large 
scale coherent field outside the disk, the poloidal field, the differential winding of this poloidal field 
leads to a large scale force-free helix that transports the magnetic energy away from the disk and 
from the dynamo. The back reaction of this force- free field (force- free except at the disk surface 
boundary) only acts as a torque on the Keplerian flow and thus the field energy of the force-free 
helix can grow at the expense of the free energy of formation of the black hole. The back reaction of 
this force-free field, being much smaller than the toroidal field, does not affect the plume formation 
by star-disk collisions. Only the much larger toroidal field affects the plumes and this in turn must 
be less than the pressure inside the disk. Thus the star disk collisions produce a robust dynamo 
where the back reaction does not quench the dynamo action at low values of field. The resulting 
exponential gain of this dynamo is an instability converting kinetic to magnetic energy. Since the 
gain is large, the dynamo fields should rapidly grow to saturation or the back reaction limit. This 
limit we conjecture is the torque corresponding to the accretion flow of angular momentum away 
from the black hole. Hence, the dynamo should convert a large fraction of the free energy of the 
black hole formation to magnetic energy. 
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A. On the Parity of Magnetic Fields 

Any arbitrary vector field C = C(r, <f>, z) can be decomposed into the sum of parts even and 
odd with respect to the reflection z — > —z, C = C e + C°. The following symmetry rules are valid 
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for an even field: 

c e r (-z) = c?(z), c;(-z) = c;(z), ct(-z) = -c e z (z), (ai) 

and for an odd field: 

C° r (-z) = -C° r (z), C° 4> (-z) = -C° 4> (z), C°(-z) = C°(z). (A2) 

Often even fields are called quadrupole type fields and odd fields are called dipole type fields. 
The last terminology reflects on the largest scale modes possible within each symmetry class and 
allows one to visualize fields of each symmetry type easily. The even and odd decomposition of an 
arbitrary field C can be performed as follows: 
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One can check that for any volume V symmetric with respect to the plane z = 

f C 2 dV= [ {C e fdV+ f (C°) 2 dV. (A4) 
Jv Jv Jv 

This implies that if C = B is a magnetic field, then the energy of the magnetic field is equal 
to the sum of the energies of its even and odd components. The even and odd components of 
solutions of equations ([3l"J) and ([35]) decouples if the mean velocity field is even, vp r {— z) = vp r (z), 
vp z {—z) = —vp z (z), z) = £l(z), the coefficient a is antisymmetric with respect to reflection 
z — * —z, and the coefficient (3 is symmetric with respect to reflection z — ► —z. Thus, even 
(quadrupole) and odd (dipole) modes will have different growth rates. The axisymmetric magnetic 
field is even if A(-z) = -A(z), B^z) = B^-z) and is odd if A(-z) = A(z), B^-z) = -B^(z). 



B. Skin Effect for the Magnetic Dynamo 

Let us consider equation (|72p written in spherical coordinates g, 6, and <f> such that 9 = 
and 9 = ir corresponds to the symmetry axis of the system. In the case of time-dependent flow 
described in section 16.11 there are no eigenmodes with a fixed frequency. Instead, the magnetic 
field can be represented as an integral over frequencies in the Fourier transformation. However, 
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in the case of a growing (and possibly oscillating) magnetic field, there is a characteristic growth 
rate T of the dynamo averaged over plume pulses. In addition, the magnetic field will possess 
oscillating Fourier components associated with the period of the emergence of plumes and, possibly, 
some intrinsic oscillatory behavior of the dynamo. We consider the behavior of one such Fourier 
component assuming the dependence A oc exp(— iut), where the complex u is the sum of the real 
and imaginary parts as uj = u' + iT. The T is the average growth rate of the dynamo, while u/ can 
take on a whole range of values, including the frequency of plumes, the Keplerian period, all its 
harmonics, etc. We impose the boundary condition for A on some sphere of radius gi n such that 
Qi n > 1 but still Qin is of the order of 1. We assume that the value of A at g = Qi n is dictated by 
the dynamo process inside Qi n . Then, for one Fourier component equation ([72]) becomes 



A 18/ 2 8A\ 1 f A 

iu} ~ = -2 7T [e IT + ~ LA > ( B1 ) 
rj g 8g \ ug J g z 



where 

L = t: ^tt. Sin 9 



sin9 89 V 99 J sin 2 9 
is the angular operator acting on A. In spherical geometry, equation (|Blj) has separable variables 
g and 9. Thus, we look for solutions in the form A = Ri(g)Qi(9) exp(— iut). 

The operator L commonly occurs in problems with axisymmetric flows, when solving the 
equation for the stream function. Since the magnetic field should be finite on the axis 9 = 0, the 
quantity 

1 d 

-(sin 9 A) 



sin 9 89 

must be finite at 9 = and at 9 = ir, because Bp = V x (Ae^). The eigenvalues and eigenfunctions 
LQi = XiQi satisfying these boundary conditions are 

A, = -1(1 + 1), Q t = sin0P/(cos0), (B2) 

where prime denotes the differentiation of the Legendre polynomial Pi(x) with respect to x and 
I = 1,2,3, .. . Besides these eigenvalues, A = is also an eigenvalue with the eigenfunction Qq = 
(1 — cos 9)/ sin 9. The first three eigenfunctions given by formula (|B2|) are 

Qi = sin9, Q2 = sin 6* cos 9, Q3 = sin# ^cos 2 9 — -J . (B3) 

The angular dependence Qi(9) determines the symmetry of the solutions. The mode proportional 
to Qo describes the radially directed magnetic field with nonzero total flux through the sphere from 
9 = to 9 = ir. All terms with / > 1 corresponds to the magnetic field with vanishing total flux 
through the sphere from 9 = to 9 = ir. The Qq term cannot be excited by the dynamo operating 
inside gi n because of V-B = condition. This is also clear from the fact that Qq — > 00 when 9 — > ir, 
which means that the vector potential cannot be well defined for a magnetic field with V ■ B / 0. 
The terms with I > 1 represent multipole expansion of the magnetic field in the far zone of the 
generation region. Ri(g)Qi(9) is a dipole term, R2(g)Q2(6) is a quadrupole term, and so on. 
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For the radial part of the solution we obtain the equation 



dg z g ag g z rj 
We introduce a new variable z = g/x where 

" CB6) 

Then, equation (|B4p reduces to the Bessel equation of imaginary argument. Solutions of this 
equation which vanishes at g — > oo are given in terms of modified Bessel function K v (z) as 



Rl = J 2z K W/z ( z ) ■ 

The Bessel functions of half -integer order can be expressed through elementary functions (e.g., 



Abramowitz Stegunl ()1972l )). Thus, we obtain for the dipole and quadrupole terms 



*M = £«-*( 1 + ;)' *m = £«-(i+§ + ? 

Finally, collecting all the terms together and retaining only the leading dipole and quadrupole 
terms, we obtain the following solution for A 

A = ai sm9^e- s/x + e~ iult + 
2g V QJ 

a 2 sin 9 cos e^e~ e/x (l + — + %-) e~ iujt , (B6) 
2g V 6 Q J 

where the coefficients a\ and a>2 should be determined by the condition of the continuity of harmonics 

of A at the surface g = gi n . The values of a\ and a 2 are determined by the dynamo action inside 

the radius gi n . We see that both dipole and quadrupole components (and all higher multipole 

components) decay as oc e~ e ' x . Using the expression (|B5p for x one obtains 

e -e/x = exp ( — |=^/v / r 2 + t,/2 + r + i^=^/v / r 2 + a;' 2 ^r) , (B7) 

where we assumed T > and oj 1 > 0. The thickness of the skin layer is determined by the real part 
of the expression under the exponent in equation (|B7j) . The larger the growth rate T, the faster the 
magnetic field decays with the radius. Also, oscillating modes with <J > decay faster with the 
radius than the steady modes with oj' = 0. Thus, far from the dynamo source one should expect 
the magnetic field to be growing in time, steadily, without oscillations. 



The characteristic length of the exponential decay of the field, l s , is found from equation ()B7p 

as 

l s = J . 2 1L= . (B8) 



For a steady magnetic field l s = y/rj/T. When g is approaching the radius of the outer boundary 
i?2) the solution (|B6P starts to "feel" the boundary condition as an ideally conducting boundary 
and the numerical results at g > R 2 are not approximated by formula (|B6|) . 
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Fig. 1. — The a — VL dynamo in a galactic black hole accretion disk. The radial component of the 
poloidal quadrupole field within the disk (A) is sheared by the differential rotation within the disk, 
developing a stronger toroidal component (B). As a star passes through the disk it heats by shock 
and by radiation a fraction of the matter of the disk, which expands vertically and lifts a fraction of 
the toroidal flux within an expanding plume (C). Due to the conservation of angular momentum, 
the expanding plume and embedded flux rotate ~ ir/2 radians before the matter in the plume and 
embedded flux falls back to the disk (D). Reconnection allows the new poloidal flux to merge with 
and augment the original poloidal flux (D). 
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Fig. 2. — Initial even (quadrupole) axisymmetric poloidal magnetic field in the simulation box. 



-52- 



x 

K) CN «- O *— CN r-o 




(b) 



Fig. 3. — The diffusion of the initial poloidal magnetic field shown in Fig. [2j The axisymmet 
poloidal field at the time t = 140 is shown. 
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Fig. 4. — Toroidal magnetic field produced by differential rotation with the Keplerian angular 
velocity VLk = r~ 3 / 2 starting from the initial quadrupole poloidal magnetic field shown in Fig. [2j 
The contours of equal magnitude of the axisymmetric toroidal magnetic field at the time t = 140 
are shown. 
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Fig. 5. — Evolution of fluxes of three components of magnetic field for the simulation presented in 
Figs. L2HH Natural logarithms of the absolute value of fluxes are plotted vs. time. Flux of B r is 
calculated through the part of the cylindrical surface r = 2 limited by lines cp = 0, 4> = n/2, z = 0, 
and z = 4. Flux of B z is calculated through the surface z = 3 limited by lines r = Ri, r = R2, 
= 0, and <p = tt. Flux of is calculated through the surface 4> = limited by lines r = R\, 
r = i?2, z = 0, and z = 4. 
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Fig. 6. — Evolution of the initial odd (dipole like) poloidal magnetic field in differentially rotating 
plasma with the Keplerian angular velocity Qk = r -3 / 2 . The poloidal axisymmetric magnetic field 
at the time t = 210 is shown. 
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Fig. 7. — Evolution of the initial odd (dipole like) poloidal magnetic field in a differentially rotating 
plasma with the Keplerian angular velocity Qk = r~ 3 / 2 . The contours of toroidal axisymmetric 
magnetic field at the time t = 210 are shown. 





Fig. 8. — Evolution of fluxes of three components of magnetic field for the simulation presented in 
Figs. [6HZ1 Natural logarithms of the absolute value of the flux are plotted versus time. The flux of 
B r is calculated through the part of the cylindrical surface r = 2 limited by lines = 0, (ft = n/2, 
z = 0, and z = 4. The flux of B z is calculated through the surface z = 3 limited by lines r = Ri, 
r = i?2, (ft = 0, and (ft = it. The flux of is calculated through the surface (ft = limited by lines 
r = Ri, v = i?2> z = 0, and z = 4. 
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Fig. 9. — Distortion of magnetic field lines of a quadrupole field by a single cylindrical plume rising 
in a fluid, which is at rest. Side view from ^-direction. The picture shows the lifting of a bundle 
of field lines. Arrows indicate the velocities of the flow at different depths in the plume. The 
corkscrew motion of the plume is clearly seen on the side view. 
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Fig. 10. — Same field configuration as in Fig. [9] viewed from the top z-direction. The picture shows 
the twisting of the bundle of field lines. Arrows indicate velocities of the flow. The time is the same 
as on Fig. [UJ 




Fig. 11. — Creation of the poloidal magnetic field from the toroidal by the rising and unwinding 
jet produced by star-disk collisions in a differentially rotating plasma. Side view from r-direction. 
Arrows indicate velocities of the flow in the frame corotating with the base of the plume. The 
picture shows the rising bundle of field lines. 
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Fig. 12. — Same field configuration as in Fig. [TT] viewed from the top z-direction. Arrows indicate 
the velocities of the flow in the frame corotating with the base of the plume. The picture shows 
the twisting of the bundle of field lines. 
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Fig. 13. — Same field configuration as in Fig. 1111 viewed from the side ^-direction. Arrows indicate 
velocities of the flow in the frame corotating with the base of the plume. The picture shows the 
formation of the loop of poloidal field lines by the rising plume. 
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Fig. 14. — Exponential growth of the dynamo magnetic field. Half of the logarithm of the B 2 
averaged over all computational domain is plotted versus time in the top- left plot. Time evolution 
of the fractions of the energy of the odd and even components of the magnetic field is in the top- 
right and bottom-left plots. The sum of the fractions of odd and even components is always equal 
to 1. 
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Fig. 15. — Time evolution of logarithms of the absolute values of the components of magnetic field 
averaged over the surfaces described in the text. Exponential growth of all three components is 
evident. 
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Fig. 16. — Time evolution of the ratio of the energy of the magnetic field in the outer domain to 
the total energy of the magnetic field in the computational domain is on the top-left plot. Time 
evolution of the fractions of energy of the odd and even components of the magnetic field in the 
outer domain is on the top-right and bottom-left plots. The sum of the fractions of odd and even 
components is always equal to 1. 
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Fig. 17. — Time evolution of logarithms of the absolute values of the radial component of the 
magnetic field averaged over the four surfaces described in the text. All four surfaces are located 
in the outer region described in the text. 
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Fig. 18. — Time evolution of logarithms of the absolute values of the toroidal component of the 
magnetic field averaged over the three surfaces described in the text. All three surfaces are located 
in the outer region described in the text. 
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Fig. 19. — Time evolution of logarithms of the absolute values of the axial component of the 
magnetic field averaged over the two surfaces described in the text. Both surfaces are located in 
the outer region described in the text. 
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Fig. 20. — Vector plots of the poloidal magnetic field of the growing dynamo at the time t = 640 
in the plane = 0. The length of arrows on the top plot is proportional to the magnitude of the 
poloidal magnetic field. Arrows on the bottom plot have unit length and are directed along the 
poloidal field. 



Fig. 21. — Three dimensional plot of the poloidal magnetic field of the growing dynamo at the 
time t = 320. The length arrows scales as the 1/3 power of the magnitude of the poloidal magnetic 
field. The magnetic field on the planes <j> = ir/2 and (j) = 3ir/2 is shown. The boundaries of the 
cylindrical region are z = ±4, R\ = 0.2, R2 = 4. 



